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1.1  Observing  Satellite  from  Aircraft 


Initiator:  J.  Buchau 

Project  No:  4643  Problem  No:  4714 

The  capability  of  program  LOKANGL  to  produce  the  position  of  a satellite  with 
respect  to  an  observation  station  was  augmented  by  adding  subroutine  CORFL. 
This  addition  allows  program  LOKANGL  to  observe  a satellite  from  an  aircraft 
as  a moving  station.  Subroutine  CORFL  calculates  for  any  moment  of  time 
the  corresponding  aircraft  position,  given  a sequence  of  its  path  segments 
and  its  flying  altitude.  Input  information  describing  the  path  of  an  aircraft 
may  be  given  either  in  the  regular  format  provided  for  stations  in  LIKANGL, 
or  in  the  format  used  by  the  researcher.  In  this  latter  case  the  special 
service  program  FLTRANS  was  created  to  transfer  such  information  to 
the  format  acceptable  by  LOKANGL.  Figure  1 reflects  the  updated  input  format. 
1.1.1  Functional  Description 

The  path  of  the  aircraft  consists  of  a number  of  segments.  Each  segment  is 
considered  a part  of  a great  circle  and  is  given  by  its  starting  and  ending 
coordinates  (longitude,  latitude),  as  well  as  the  flight  altitude.  Each 
flight  segment  is  thus  given  as  two  stations.  Subroutine  CORFL 
interpolates  the  aircraft  position  at  any  desired  instant,  and  this  location 
is  then  considered  by  LOKANGL  as  an  observation  station  position  for  output 
calculations.  When  a number  of  flight  segments  is  given,  program  LOKANGL 
should  be  provided  with  a corresponding  number  of  pairs  of  station  cards  for 
starting  and  ending  positions.  Number  of  stations  in  LOKANGL  input  should  be 
equal  to  the  number  of  flight  segments. 


BEST  AVAIWBIE  COPY 


INPUT  F3PN»I  FO<  PPOliMM  LOKANCL 


0»T»  OEC“’  SilUP 

CAPO  -aLS  JtSCPIPTION 

_1_ _1 COPE  OF  orbital  DFTERHINATION  FQPH  

l = FORH  NO.  } » SCF  2-CARO  POS.-VFl."  ^T" 

?*FONH  NO.  ? I AOC  ?-CAPO  ELENFNT  PATA  SET 

3=3-CAR0  ELENENT  OATA  SET 

4=osculatino  elements 

S:FORN  NO.lt  AOC  S-CARO  OATA  SET 

ONE^b^NORE  ELEMENT  S^fS  * " 

f-1  X CR  U ON  FIRST  CARD  CF  ALL  SU8SF0IIFMT  FLEMENT  SETS 

1 1-1  NO.  OF  STATIONS.  IF  NO  STATIONS  USF  D 

f CODE  I OR  0 FOR  PRINT  CONTROL  OF  STATION  OATA 

_0  = PRI^NT  BY  TI^  ONLY 

1 = PRINT  BY  SfATIONS' 

T-io  o = stanoaro 

IONOHT  = SU8-IONOSPHERIC  altitude  IfHI  13 
13-lS  MINELV=  NIN  ELEVATION  FOR  SAT  VIFhINF.  T3 

T» STATION  LOCATION  CMOS, IF  NMS.GT.O 

I- S  I.'C.  OF  STATION  «NUN3FPI 

T-8  ?=GE00eTir  SYSTEN 

-I=INOICATES  MOVING  STATION  (AIorpAril 
=-?3  STATION  GEODETIC  LATITUDE.  DEGREES 
?4-3J  STATION  LONGITUDE  (POSITIVE  HESTI . OEGoFES 

STATION  HEIGHT.  METERS  _ 

fl-T?  NAME  OF  StATION 

4 1-15  TIME  INCREMENT  IN  SECCNDS.  F15.5 

lT-i9  TIME  interval  FOP  PRINT  OP  TAPE  ^ VPITF 

START  COLSI  IT- 18 1 30- 21 1 23-241 26- 79:  31- 34;  36-41 
_FIN^  COLSI  41-44;46-47;49-5aL52-‘=S;'7-60;62-6T 
TIME  month;  "tay ; year;  hcio;  mtn;  sec 

FORMAT  iz;  12;  12;  F3.  i;  FT.i;  F5.3 

9 CODE  FOR  print-out;  0=NO  print;  t=POIMT  FOLLOHINGt 

7 1=POSITION  ANU  VELOCITY 

4 l=SUa-SATELLITE  DAT! 

9 1 = 0SCULATING  ELEMENTS 

• 1 = 09SERVATI0N  OATA 

p-10  0=STANDARO  BINARY  ON  TAPE3 

-l=°ERIGEE-APn&EE  DATA  ON  TAPE/ 

II- 13  l.O-SINGLE  PASS  PEFRAC.  CORRECTION  tSTANOARPI 

8 ...t-I'!.  “END  0F_PRp8LEM“ 

A lOVE  MAY  3£  PEPEATEO  N TIMES 
LAST  CARTi  "9"  IN  COL.  1 

COM«rNT  II  USAGE  OF  MOVING  OOSFRVATION  STATION  SthulATING 

AIRCTACT  REQUIRES  INPUT  OF  SEVERAL  (N.ZI  CON'EEUTTVE  AIRCRAFT 
PDSI'ION  rOOPDINATE  CAROS 

riMMFNT  7t  EACH  FLIGHT  SEGMENT  MIST  PE  SUPPLIED  MITH  T MO 
POSITION  COORDINATE  CARDS  SIGNIFYING  STaRT  AND  ENr  OF 
SEGMFNT.  for  each  pair,  NUMOER  of  STATIONStNMSt  TS 
INCPrisro  HY  ONE  I NMS  = NMS A1 1 
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Figure  1.  Input  Format  for  Program  LOKANGL 
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Input  format  provided  by  the  initiator  consists  of  a number  of  cards,  each 
of  which  provides  the  coordinates  and  times  of  turning  points  along  the 
flight  path.  Program  FLTRANS  uses  this  information  to  create  starting  and 
ending  points  for  each  flight  segment,  as  well  as  all  other  parameters 
required  by  LOKANGL. 

1.1.2  Mathematical  and  Logical  Procedures 

Given  the  total  time  of  flight  T along  the  flight  segment,  starting  and 
ending  coordinates  ((1)2X2)  and  x^)  of  the  flight  segment,  height  of  the 
flight  H,  and  earth  radius  R;  the  position  of  aircraft  (()),X)  at  the  time  t is 
determined  as  follows: 


1.  X = arc  Cos(Sin  ({)2Sin  Cos  (t)2Cos  (t)2Cos(X2-X2) ) . 

2.  A = arc  Cos{(Sin  (t)2-Cos  X Sin  ())2)/(Sin  X Cos  (1)2)). 

3.  X = Xt/T. 

4.  (j)  = arc  Sin  (Sin  (j),Cos  X + Cos  ())j^Sin  xCosA) 

5.  X = Xi+  arc  Sin(Sin  A Sin  x/Cos  (()). 


I 
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1.2  Long  Term  Ephemeris  for  Near-geostationary  Satellites 

Initiator:  J.  Klobuchar 

Project  No:  4643  Problem  No:  4887 

The  orbits  of  near-geostationary  satellites  are  significantly  influenced 
by  resonant  perturbations  due  to  the  lower  order  tesseral  terms  in  the 
geopotential,  and  to  a lesser  extent  also  by  the  luni-solar  gravitational 
forces.  These  effects  show  up  only  over  a period  of  months,  and  some 
analyses  are  available  in  the  literature  for  the  lowest  tesserals 

(J2  2 3nd  33^3)  and  for  low  eccentricity  and  inclination.  Luni-solar  forces 
are  of  course  time  dependent.  The  objective  of  this  work  was  to  improve  the 
capability  of  the  AFGL  satellite  orbit  generation  programs  to  predict  long 
tif?rm  longitudinal  drift  for  geostationary  satellites. 

A LOKANGL-type  program  which  generates  orbits  based  on  elements  and  their 
rates  of  change  is  suitable  only  when  pronounced  effects  are  involved  such 
as  those  due  to  the  J2  and  J3  zonal s for  low  altitude  polar  orbits.  On  the 
other  hand,  the  elaborate  analysis  required  to  calculate  geostationary 
resonances  is  implicit  in  program  ROPP  which  includes  tesserals  thru  c. 

b >b 

lunar  and  solar  perturbations,  and  an  atmospheric  drag  model. 

Program  ROPLOK  was  accordingly  developed  by  removing  the  drag  calculations 
of  ROPP,  simplifying  to  a no-overlay  program  capable  of  executing  in  77K 
(octal),  and  with  input  similar  to  LOKANGL.  (See  Figure  2.  LOKANGL-type 
output  not  available).  Accurate  knowledge  of  the  mean  motion  and  the  semi- 
major axis  is  essential  to  predict  the  long  term  resonance  amplitude  in 
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nopir  pmeoiction  pro&ram  roplok 


0»TA  nF'':<  SITUP 

NAHFLISr  INPUT  FOR  CHECK-OUT  HAY  PRECEDE  NO»HAL  TNPUT 
THI*;  JATA  IS  ENTERED  HETHEEN  tlNPUT  AND  »YNn  CAROS 

FARO  -DLi  OESCRI»TION 

1 1 orpital  oeterhinatiom  form  code 

INFORM  NO.  S I SCF  2-CARD  POS.-WFL.  FFT 

Z=FORM  NO.  2 t A3C  2-CARD  ELEMENT  DATA  SET 

5 = FORM  NO,  l_l MC  5-CARO  DATA  «> '"T 

ONE  OR  TMO  ELEHCNT  SETS 

A-T  X OR  U ON  first  CART  CF  EACH  ELEMENT  SrT 
NOTFJ  IF  TMO  ELtHENT  SETS  ARE  INP'J», 

MEAN  NOTION  OF  THE  FIRST  SET  IS  AnjIlS^'P 
TO  MATCH  TRUE  ANOMALY  OF  THE  SECON P SET. 

T 1-T  NO.  OF  STATIONS.  IF  NO  STATIONS  USE  0 

F CODE  1 OR  0 FOR  PRINT  CONTROL  OF  STATTON  DATA 

0 = PRINT  3Y  TIME  OriLY 

1 = PRINT  3Y  STATIONS 
’llD  O^STANOWO 

lONOHT  =SUB-I0N0SPMERIC  ALTITUDE  tKM)  13 

5*  STATION  LOCATIOr  CARQS.IF  NNS.GT.O 

I- S  I.O.  OF  STATION  INUM3ER1 

II- 25  STATION  OEOOEI It  LATITUDE.  OEGREFS 

^^-'•9  STATION  LONGITUDE  (“JSITIVE  EASTT.  DFE.PFES 
<•1-55  STATION  MEIGHT.  METERS 
fl-22  NAME  OF  STATION 

<•  1-5  TIME  INCREMENT  IN  'tliUlCS.  15 

11-42  TIME  span  for  point  OF  TAPE3  MPITE 

start  COLSI  11-12;14-15;1T-19;20-?11?3-?4;2P-?T 

final  CX51  31-32;  T4-3E;3T-38;40-41!fc3-44;4fe-4T 
TIME  YEARtMONTH;  DAY!  HOIlo;  MTN;  SEC 

FORMAT  12;  12;  12;  I’:  T2;  12 

CODE  FOR  print-out;  0 = NO  print;  1=poT*IT  FOLLOMING; 
-..IfPQSLLlON  A_NO  VFLOCITY 
4 l = SUB-SATELLlTt  DATA 

F i=mean  keplerian  elements 

F IsCPSERWATION  DATA 

?-ia  TAPE  3 control;  stanoaro=o 

a IjiJJi.  -•:£NQ_flELPH05tEM~  

A OQV;  MA>'  3£  ‘CREATE  3 N TIMES 
LAST  SI’-:  -i"  IN  SOL.  1 


Figure  2.  Input  Format  for  Program  ROPLOK 


BEST  AVAIWBLE  COPY 


longitude.  This  facility  was  achieved  by  accepting  and  adjusting  to  two 
sets  of  elements  which  would  ideally  be  provided  separated  by  a number  of 
weeks.  Subsequently,  the  tesseral  selecting  subroutine  RESON  was  also 
incorporated  into  ROPP. 

Subroutine  RESON  determines  whether  the  net  motion  specified  by  the  elements 
approaches  any  resonance  commensurability  IB/IA  revolutions  per  day,  where 
lA.IB  may  be  any  integers  from  1 thru  6.  If  such  a ratio  of  integers  is 
found,  a search  is  conducted  for  all  L,M  tesseral s where  the  allowable  P,Q 
subscripts  result  in  significant  geopotential  contributions.  Thus,  solve 
for  L,M,P,Q  from 

IB(L-2P+Q)  = IA(M) 
for  the  following  ranges; 

L=2  thru  6 
M=0  thru  L 
P=0  thru  L 

Q=-l  thru  1 for  eccentricity < .08 

-2  thru  2 for  eccentricity < .16 

-3  thru  3 for  ecentricity  ^.16 

The  corresponding  inclination  function  FLMP  and  eccentricity  function  GLPQ 
are  obtained  using  existing  routines  FINK  and  FECC  respectively. 

The  contributing  geopotential  field  is  then  given  by 

VLMPQ  = — ^ •J,  ^-FLMP-GLPQ 

where  u = 398603.15  km^/sec^ 

and  a is  the  semi-major  axis  in  earth  radii. 
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Program  ROPP  is  set  to  always  include  the  J2  and  zonal s because  of  their 
usual  importance.  In  the  sample  tabulation  below  VLMPQ  values  are  also 

Q 

listed  for  all  cases  greater  than  10  for  a satellite  with  mean  motion 
n=l. 00273  (and  a=6.6109),  e=0. 00907,  i=2.50. 


L 

M 

P 

Q 

VLMPQ 

* 

2 

2 

0 

0 

+7.34x10"^ 

★ 

3 

1 

1 

0 

-6.09x10"^ 

3 

2 

0 

-1 

-2.32x10"^ 

3 

2 

1 

1 

-6.95x10"^ 

* 

3 

3 

0 

0 

+6.75x10"® 

★ 

4 

0 

2 

0 

+1.93x10"® 

4 

1 

1 

-1 

-1.56x10"® 

4 

1 

2 

1 

+7.81x10"® 

★ 

4 

2 

1 

0 

-3.44x10"® 

4 

3 

0 

-1 

-5.61x10“® 

4 

3 

1 

1 

-1.68x10"^ 

★ 

4 

4 

0 

0 

+2.21x10"® 

5 

1 

2 

0 

+3.52x10“^ 

5 

2 

2 

1 

+1.77x10"® 

* 

5 

3 

1 

0 

-2.70x10"® 

* 

5 

5 

0 

0 

+5.37x10"® 

6 

0 

3 

0 

+1.43x10"^ 

6 

2 

2 

0 

+3.52x10"^ 

For  a selection  level  of  1.0x10"®  the  terms  starred  are  significant;  a sample 
of  the  corresponding  print-out  is  shown  in  Figure  3. 
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Figure  3.  Selection  of  Resonance  Terms  in  Program  ROPLOK 


r " 

In  order  to  gauge  the  relative  long  term  effect  of  these  terms,  a number  of 
different  combinations  were  run  over  a 200  day  span  for  the  above  satellite, 
and  for  a satellite  with  e=0,01  and  i=5°.  Table  1 shows  longitude  drift 
predictions.  With  no  extra  terms,  the  drift  is  small.  As  expected,  inclusion 
of  2 shows  the  resonance  drift  well,  with  the  ^ and  ^ tesserals 
essentially  completing  the  picture.  The  remaining  terms  contribute  little. 

The  selection  threshold  for  VLMPQ  is  adjustable  in  ROPP,  but  the  default  value 
is  set  at  1.0x10"^. 

Finally,  a study  was  made  for  future  use  on  a SCATHA  type  satellite.  For 
n=l. 018965,  e=0. 17907  and  i=2.5°  the  satellite  drifts  eastward  at  6°/day.  Thus 
no  resonance  commensurability  exists,  but  the  same  combinations  of  terms  as 
before  were  forced  in  for  testing.  After  200  days  a minor  effect  of  including 
the  (2, 2, 0,0)  term  was  evident,  giving  a longitude  1.5°  East  of  the  longitude 
with  no  extra  terms.  Additional  terms  produced  no  tangible  effects.  Since 
the  higher  mean  motion  tends  towards  the  (2, 2, 0,-1)  term,  this  was  also  tested. 
The  CP  execution  time  was  unexpectedly  quadrupled,  but  again  there  was  no 
effect  on  the  longitudinal  drift.  A SCATHA  type  satellite  is  not  near- 
geostationary, and  additional  terms  are  properly  excluded. 
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GEOSYNCHRONOUS  SATELLITES 


L.M.P.Q  TEST  (200  day  run) 


Terms 

Included 

(2,0), (3.0) 
always  present 

Initial 

All  (2,2)(3,1)(3.3)(4,0) 

(4.2) (4,4)(5.3)(5,5) 

(2. 2) , (3,1), (3, 3) 

(2.2)  alone 

(2. 2) . (3.1) 

(2. 2) . (3, 3) 

No  Extra 


n=l. 002730 

i=2.5,e=. 00907  i=5.0,e=.01 
East  Longitude 


285.86 

285.81 

266.51 

266.63 
(also  3,2) 

266.65 

266.76 

262.10 

262.25 

263.29 

263.41 

265.43 

265.56 

288.35 

288.43 

Table  1.  Longitude  Drift  Due  To 
Selected  Geopotential  Terms 
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1.3  Subroutine  ECLIPS 

Initiator:  E.  Robinson 

Project  No:  0001  Problem  No:  4517 

Subroutine  ECLIPS  was  written  to  be  called  by  rapid  orbit  prediction  program 
ROPP  to  produce  in-out  eclipse  times  of  a satellite.  Subroutine  ECLIPS  was 
developed  in  order  to  replace  the  existing  subroutine  which  produced  in-out 
eclipse  times  using  iteration  methods  which  required  unnecessarily  long 
execution  times  and  was  not  able  to  produce  smooth  and  reasonably  exact 
in-out  eclipse  time  calculations.  The  new  subroutine  ECLIPS  is  based  on 
an  analytical  approach  developed  in  chapter  5 of  Reference  7 which  resolves 
the  above  difficulties. 

1.3.1  Functional  Description 

Subroutine  ECLIPS  determines  the  entrance  and  exit  true  anomalies  of  the 
satellite  into  the  shadow  of  the  earth,  and  the  duration  of  eclipse.  It 
assumes  no  flattening  of  the  earth  and  no  shift  of  the  earth  in  its  orbit. 

Umbra  and  penumbra  effects  are  also  neglected.  These  assumptions  will 
result  in  a quartic  equation  in  the  cosine  of  the  entrance  and  exit  true  anomalies. 

A given  orbital  ellipse  is  fully  defined  as  to  dimensional  size  and  orienta- 
tion in  the  equatorial  or  right  ascension-declination  coordinate  system  by 
semi-major  axis-a,  eccentricity-e,  incl ination-i , longitude  of  the  ascending 
node-J<,  argument  of  perigee  oj  and  the  time  of  perifocal  passage  T.  Given 
the  above  information,  the  shadow  function  ^ is  calculated  and  composed  into 
the  quartic  equation  for  determination  of  entrance  and  exit  true  anomalies 
of  the  satellite  into  the  shadow  of  the  earth.  To  solve  this  equation 
subroutine  POLRT  is  used. 
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Based  on  true  anomalies,  the  mean  anomalies  for  entrance  and  exit  from  the 
shadow  and  therefore  in-out  times  are  calculated. 


1.3.2  Mathematical  and  Logical  Procedures 


Given:  a,e,i,fi,HT  we  calculate  first  the  orthogonal  set  of  vectors 


Px  = Cos  03  Cos  -Sin  0)  Sin  Cosi 


Py  = Cos  (i)  Si  n +Si  n w Cos  SI  Cos  i 


(£  is  taken  as  pointing 
toward  peri  focus.) 


P^  = Sin  u)  Sini 


= Sin  0)  Cos  SI  -Cos  o)  Sin  SI  Cosi 


Qy  = Sin  u)  Sin  SI  +Cos  03  Cos  SI  Cosi 


= Cos  03  Sin  i 


(2  is  in  the  orbit  plane 
advanced  to  P by  a right 
angle  in  the  direction  of 
increasing  time  anomaly.) 


We  calculate  then 


X P +Y  P +Z^P 
sx  sy  sz 

p 2 2-“— z — r 

(X  +Y  +Z 


X^q^+YgQy+ZgQ^ 
C = , 7—^ 2 57 


(X  3.Y  ^.Z  3)' 


where  X , Y , Z are  coordinates  of  the  sun. 
s s s 


Shadow  function  of  the  Earth  is  produced  by 


2 22  22 

S = a (1  + e Cos  0 ) + P (3Cose  + ?Sin0  ) - P 


where  S = 0 is  the  condition  for  entrance  or  exit  from  the  shadow  and 


P = a(l-e  ) is  the  semi -parameter  of  the  orbit.  0 is  the  true  anomaly. 


Satellite  enters  the  shadow  if  S changes  sign  from  to  "+".  Exit 
from  shadow  is  determined  by  changing  of  sign  of  S from  "+"  to 
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Rejection  of  spurious  roots  is  accomplished  by 
Cos  <1^  = BCose  + 5 Sin  0 < 0 

since  il>  is  the  angle  between  the  radius  vector  to  the  sun  and  the  radius 
vector  to  the  satellite  at  entrance  or  at  exit  points. 


Solution  of  the  equation  S = 0 

A Cos'*©  + A Cos^e  + A Cos^e  + A Cose  + A =0 

0 12  3 4 

with 

A^  = [(ag/p)-e‘*-2(ag/p)(^'-e')e2  + (6^  + 

A^  = j^4(ag/p)‘'e^  - 4(ag/p)^  (^^-6^)ej 

A^  = [^6(ag/p)‘*e"  - 2{ag/p)"  - 2(ag/p)'(  1-C')e=^+ 

2(£"-B")(l  -C")  - 46^5^] 

A^  = 4(ag/p)‘*e  - 4(ag/p)^(l  -5^)ej 

A,  = fcte/P)'*  - 2 ^ 


25 


REFERENCES 

1)  Wexler,  D.M.,  "Rapid  Orbit  Prediction  Program  (ROPP)", 

TRW-09967-6001 -RO-00,  prepared  under  AFCRL 
Contract  No.  F19628-68-C-01 70,  October  1968. 

2)  Sparks,  S.T.,  "Users'  Manual  for  the  TRW  Rapid  Orbit  Prediction 

Program",  TRW-09967-6002-R0-00,  prepared  under 
AFCRL  Contract  No.  F19628-68-C-01 70,  October,  1968. 

3)  Kaula,  W.M.,  Theory  of  Satellite  Geodesy,  Blaisdell  Publishing 

Co.  (1966) 

4)  Blitzer,  L.;  Boughton,  E.M.;  Kang,  G.;  and  Page,  R.M., 

"Effect  of  Ellipticity  of  the  Equator  on  24-Hour 
Nearly  Circular  Satellite  Orbits",  J.  Geophysical 
Research,  Vol . 67,  No.  1,  January  1962. 

5)  Allan,  R.R.,  "Perturbations  of  a Geostationary  Satellite  by  the 

Longitude-Dependent  Terms  in  the  Earth's  Gravita- 
tional Field",  Planetary  Space  Science,  Vol.  11, 
1963. 

6)  Allan,  R.R.,  "Satellite  Resonance  with  Longitude-Dependent 

Gravity-II",  Planetary  Space  Science,  Vol.  15, 

1967. 

7)  Escobal,  p p ^Methods  of  Orbit  Determination,  John  Wiley  & Sons, 

Inc.,  N.Y ./London/Sydney  (1965) 


26 


Section  2,  Ionospheric  Research  Support 


Author:  B.  J.  Guz 
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Initiator:  J.  Buchau 


Project  No:  4643  Problem  No:  4714 

Program  WRLDMAP  is  designed  to  provide  the  user  with  12  or  30-inch  wide  plots 
of  the  land  mass  outlines  within  desired  limits.  These  limits  are  represented 
by  geographic  or  geomagnetic  longitude  and  latitude.  Together  with  a reference 
longitude,  this  provides  full  information  for  plotting  a map  of  specific  areas 
of  either  hemisphere  of  the  Earth.  Program  WRLDMAP  is  able  to  plot  maps  with 
a given  scale  up  to  the  maximum  size  allowed  by  the  plotting  facilities.  The 
geographic  or  corrected  geomagnetic  coordinate  backgrounds  can  be  superimposed 
if  desired.  A sample  plot  is  reproduced  in  Figure  1. 

2.1  Functional  Description 

Input  for  the  program  consists  of  particular  values  of  longitude  and  latitude, 
factors  indicating  the  scale  of  maps,  and  the  desired  coordinate  background, 
either  geographic  or  geomagnetic.  Table  1 shows  the  card  input  format. 

After  initialization,  the  continental  outline  data  file  is  used  and  the  land 
mass  outlines  are  plotted  within  specified  bounds.  Then  subroutine  SUBMAP 
plots,  if  desired,  the  coordinate  lines  of  geographic  or  geomagnetic  longitude 
and  latitude.  The  plot  is  supplied  with  a framed  label. 

Program  WRLDMAP  possesses  four  options  for  confining  the  geographic  or  the  geo- 
magnetic plot  with  the  geographic  or  the  geomagnetic  coordinate  background. 

The  program  has  the  capability  of  plotting  a map  for  each  data  card  provided. 
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A blank  data  card  terminates  the  job.  This  program  requires  magnetic  tape  or 
disk  files  of  digitized  land  mass  outlines,  and  the  Hakura  or  Gustafsson  table 
of  geomagnetic  coordinate  conversions. 

2.2  Mathematical  or  Logical  Procedures 

The  proper  data  (permanent)  files  to  be  read  into  the  array  AAA  are  determined 
by  the  values  of  K3.  FACTOR  determines  the  proper  plot  size,  WRLD  determines 
whether  the  land  mass  outlines  shall  be  plotted.  If  the  land  mass  coordinate 
(X,  Y)  is  such  that  X or  Y is  greater  than  900,  the  coordinate  is  ignored. 

(Tape  2 and  Tape  3 coordinates  were  obtained  from  the  Geophysical  Institute  of 
the  University  of  Alaska  and  are  compatible  with  their  use  of  999.9998  as  a 
"non-coordinate"  abcissa  or  ordinate).  CORRGM  is  called  if  and  only  if  K3  is 
equal  to  1 in  order  to  convert  from  geographic  to  corrected  geomagnetic  coor- 
dinates when  necessary. 

If  (X,  Y)  is  the  "ordered  set"  of  coordinates  under  consideration  for  plotting, 
the  plotting  pen  is  lifted  if  and  only  if  SQRT((X.-X^_j)^  + (Y.-Y.  j)^)>.l,  or 
(X^._P^  + (Y._j^)^>100;  and  the  pen  is  moved  if  and  only  if  X.^+V(^00.,  where: 

X.  = X|^  cos  y,^,  Yj  = Xj^siny^j  , X^^  = (90.-Xg)/5., 

Y.  = AMOD  (y,  - 90.,  360)*  TT/lSO. , and  (x,,  y,)  is  the  geographic  or 

Da  ad 

corrected  geomagnetic  coordinate  (in  degrees)  obtained  from  Tape  2 (or  Tape  3) 
or  CORRGM. 

The  pen  moves  only  to  (-X^. , Y^. ) for  a north  or  south  (respectively)  polar  view 
of  the  hemisphere  determined  by  the  values  of  K1  and  K2. 

In  subroutine  SUBMAP,  the  geographic  or  corrected  geomagnetic  coordinate  system 
map  determined  by  the  values  of  Kl,  K2  and  K3  is  then  superimposed  on  the  exist- 
ing polar  coordinate  system  if  the  value  of  MAP  is  not  "NO". 
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Card 

Variable 

Card 

No. 

Name 

Col. 

Format 

Variable  Description 

1 

K1 

1-3 

13 

"1"  for  northern  hemisphere 
"40"for  southern  hemisphere 

K2 

4-6 

13 

"50"for  northern  hemisphere 
"89"for  southern  hemisphere 

K3 

7-9 

13 

”1"  for  geomagnetic  bounds  on  geographic 
background 

"2"  for  geographic  bounds  on  geographic 
background 


VIEW 

11 

A1 

WRLD 

13-14 

A2 

MAP 

16-17 

A2 

BEAT 

21-17 

F7.2 

ELAT 

28-35 

F7.2 

BLON 

36-43 

F7-2 

ELON 

44-51 

F7.2 

CLON 

52-53 

F7-2 

FACTOR 

60-67 

F7.2 

FACT 

68-75 

F7*2 

"3"  for  geomagnetic  bounds  on  geographic 
background 

"4"  for  geographic  bounds  on  geographic 
background 

"N"  for  north  polar  view 

"S"  for  south  polar  view 

"NO"  to  suppress  plot  of  land  mass  outlines 

"NO"  to  suppress  plot  of  coordinate 
background 

Beginning  latitude  of  plotted  region 
Ending  latitude  of  plotted  region 
Beginning  longitude  (east)  of  region 
End  longitude  (east)  of  the  region 
Reference  longitude-vertical  on  the  plot 
"1"  for  12  inch  plot 
"2"  for  30  inch  plot 

"0"  indicates  that  region  will  be  plotted 
in  the  maximal  possible  scale,  otherv/ise 
given  scale  should  be  entered  as  degrees 
latitude/inch. 


( 


Table  1.  Card  Input  Format  for  Program  WRLDMAP. 


I 

i 


i 
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The  major  part  of  subroutine  SUBMAP  is  devoted  to  calculating  the  intersection 
points  of  the  boundaries  of  the  plotted  region  with  lines  of  constant  longitude 
and  latitude.  Since  these  liries  could  intersect  with  more  than  one  boundary 
line  and  come  from  unpredictable  directions,  fairly  complicated  logic  is  involved. 
Algorithms  for  calculating  polar  coordinates  for  intersection  points  consists  of 

four  parts.  They  check  whether  the  intersection  of  the  given  longitude  (lati- 
tude) line  with  one  of  four  bounds  ("upper"  and  "lower"  latitude  bound  and 
"right"  and  "left"  longitude  bound)  does  exist.  For  this  purpose  a value  is 
calculated  showing  whether  the  two  consecutive  points  of  latitude  (longitude) 
line  are  on  both  sides  of  the  boundary  or  not.  Then,  if  intersection  with  a 
"latitude"  bound  has  occurred,  the  intersection  points  (polar  coordinates)  are 
calculated  by  the  formulas 

P -P  P-Pi 

6 = 0,-^ + 8o  and  P = RMAX  or  P = RMIN 


where  RMIN  and  RMAX  are  beginning  and  ending  boundary  latitude.  For  the  in- 
tersection with  "longitude"  boundaries  the  following  formulas  are  used. 

Sin  (02  - 0^ 


and  0 = 0MIN  or  0MAX 


Pj  Sin  (0-0j)  + P2  SIN  (©2-0) 


where  0MIN  and  0MAX  are  beginning  and  ending  boundary  longitude. 
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3.  Geopotential  Model  Term  Selection  For  Satellite  Orbit 
Computations 

Initiator:  K.  S.  W.  Champion 

Project  No:  427M  Problem  No:  4867 

3.0  Introduction 

Among  the  most  time  consuming  tasks  in  most  modern  earth  satellite  orbit 
computations  is  the  computation  of  the  earth's  gravitational  attraction. 

This  is  because  geopotential  models,  as  derived  from  analysis  of  artificial 
satellite  tracking  data,  surface  gravity  data  and  other  data,  contain  up  to 
several  hundred  terms  (usually  spherical  harmonics).  An  analysis  of  the 
perturbing  effect  of  these  terms  on  a satellite  orbit,  such  as  that  carried 
out  by  Kaula^  shows  that  for  any  given  satellite  only  a few  terms  will  be 
important.  Furthermore  atmospheric  drag  (for  low-altitude  satellites)  and 
in  some  cases  tracking  data  errors  (noise,  biases  such  as  station  location 
errors,  etc.)  will  overshadow  the  effects  of  many  geopotential  terms^’^. 

Thus  it  seems  worthwhile  to  develop  a selection  procedure  for  particular 
satellites  to  select  only  the  most  important  terms  of  a given  geopotential 
such  that  the  remaining  terms  will  not  have  an  appreciable  effect  on  the 
accuracy  in  orbit  determination  and  prediction. 

In  this  report  we  therefore  give  an  error  analysis  of  the  effect  of  model 
changes  on  determination  of  orbit  parameters  from  tracking  data  and  the 
subsequent  effects  on  predicted  orbit  positions.  This  is  used  with  various 
simplifications  to  select  a subset  of  terms  from  a given  geopotential  model 
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to  yield  a specified  r.m.s.  predicted  position  difference  relative  to  the  given 
model  at  a specified  time  beyond  the  tracking  fit  span.  Applications  are 
given  for  3 low-altitude  satellites  and  one  medium  altitude  satellite  using 
the  Smithsonian  Standard  Earth  III  (SEIII)  geopotential  model. ^ 

Slowey,  et  al , have  also  given  a selection  procedure  which  uses  only  the 
first-order  perturbation  amplitudes  on  the  orbit  due  to  individual  harmonics 
without  relating  these  to  prediction  accuracy. 

3.1  Analysis 

The  analysis  of  the  effects  of  individual  harmonic  terms  on  satellite  orbit 

3 

determination  is  given  elsewhere.  Briefly,  each  spherical  harmonic  gives 
rise  to  a series  of  secular  and  oscillating  perturbations  to  a given  un- 
perturbed orbit. ^ The  secular  contributions  are  limited  to  zonal  harmonics 
(those  independent  of  longitude).  The  frequencies  and  amplitudes  of  the 
periodic  contributions  can  be  computed  in  straight-forward  fashion.  Although 
the  number  of  such  terms  is  infinite,  one  may  obtain  a finite  series  by 
limiting  the  expansion  to  low  powers  in  the  eccentricity. 

The  effect  of  each  perturbation  on  a satellite  orbit  obtained  by  least 
squares  fitting  of  tracking  data  are; 

1)  The  change  in  an  orbit  of  specified  initial  conditions, 
which  can  be  computed  using  the  theory  of  Reference  1. 

2)  The  change  in  the  orbit's  initial  conditions,  which  can  be 
computed  from  a least-squares  fit  to  the  perturbation. 

Considering  only  the  in-track  component  of  the  orbit,  assuming  a second 
degree  polynomial 
3 

s(t)=E  P.-t^ 
i = l 


r 


where  represents  the  original  orbital  position,  the  original  mean  motion. 


and  P-  the  time  derivative  of  mean  motion  (as  produced  by  drag),  and  assuming 


uniform  time  distribution  of  the  tracking  data,  we  obtain,  for  the  changes 
dl  due  to  a periodic  orbital  perturbation  A cos  (at  + q^): 

1^1 2 (2  sin 


dp  ^ 

aT 


in  q -3  sin  Q ) i (2  cos  q + 3 cos  q ) 

aT  ® 


dP^  = 


aT' 


[ 


■60 


wry 


(sin  q - sin  q ) 


2 I 12  (2  sin  q - 3 sin  q^ ) (168  cos  q + 192  cos  q ) 


aT 


dP^  »^_A_ 


aT- 


[ 


-360  ( sin  q - sin  q ) 

o • 0 

(aT)^ 


t 


30(  sin  q - sin  q ) + 180  ( cos  q + cos  q ) 
° aT  ° 


-360 


riTf)' 


(sin  q - sin  q ) 
o 


]. 


where 


T = length  of  fit  span, 


q = aT  + q. 


From  this  the  r.m.s.  prediction  error  at  any  time  can  be  plotted  as  a function 
of  the  ratio  of  the  fit  span  to  the  perturbation  period,  A series  of  maxima 
are  found,  the  largest  occurring  for  fit  span  nearly  equal  to  perturbation 
period.  In  particular  very  long  period  resonances  do  not  always  produce  the 
largest  effects. 

3.2  Applications 

The  preceding  results  were  applied  to  the  SEIII  geopotential  for  satellites 
DB-7,  DB-8,  DB-9  and  OVl-17  (See  Table  1).  The  longer  fit  span  and 
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prediction  time  used  for  OVl-17  is  customary  for  the  higher  altitude 
satellites  with  less  drag. 

Radial  and  cross-track  amplitudes  due  to  perturbations  in  semi-major  axis, 
eccentricity  and  inclination^  were  added  vectorially  to  the  in-track 
components  to  produce  the  total  error  estimates  shown  in  Tables  2 and  3. 

For  low  altitude  satellites  it  was  found  that  periodic  perturbations 
in  the  perigee  height  of  the  density  scale  height)  due  to  perturbations 
in  the  semi -major  axis,  argument  of  perigee,  and  eccentricity  could  sometimes 
produce  significant  effects  on  in-track  error,  so  these  have  been  included. 

The  geopotential  model  was  reduced  in  each  case  by  rejecting  the  terms  with 
the  smallest  errors  such  that  the  sum  of  squares  of  their  errors  is  less 
than  the  square  of  the  indicated  criterion.  Thus  we  aim  to  limit  the  r.m.s. 
error  at  a specified  prediction  time  to  a specified  tolerance.  The  terms 

t 

; remaining  (those  not  rejected)  are  underlined  in  Tables  2 and  3.  Computations 

I were  carried  out  with  program  TESAMP,  originally  written  to  compute 

i 

I first-order  tesseral  perturbation  amplitudes.^  This  program  was 

modified  to  include  the  computations  described  here.  A total  of  12  cases 

; for  the  3 DB  satellites  was  processed  using  program  CELEST^  and  Doppler 

Beacon  data.  CELEST's  geopotential  routine  was  modified  to  take  better 
account  of  the  structure  of  the  selected  geopotentials,  as  in  its  original 
I form  termination  of  recursive  computational  loops  is  determined  solely  by 

the  maximum  L index  and  no  tests  are  made  for  vanishing  coefficients.  Since 
maximum  L is  nearly  as  large  for  the  selected  as  for  the  original  geopotential 
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Table  2.  R.M.S.  Prediction  Error  Effects  in  Meters  for  Terms  of  SEIII, 

as  Function  of  Their  Subscripts  L,  M;  Satellites  DB-7,  DB-8,  DB-9 
(Selected  Terms  Underlined) 
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Table  3.  R.M.S.  Prediction  Error  Effects  in  Meters  for  SEIII,  • 
Satellite  OVl-17  (Selected  Terms  Underlined) 
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little  savings  is  gained  in  CP  time.  Improvements  of  over  50%  were  gained 
however  by  early  termination  of  recursive  loops  and  avoiding  unnecessary 
computation  for  zero  coefficients. 


OVl-17  skin-track  data  from  ADCOM  was  used  with  program  CADNIP  for  fits, 
and  BADMEP^  for  prediction  errors.  Unlike  CELEST  these  programs  test  for 
zero  coefficients  and  hence  the  geopotential  computation  time  is  proportional 
to  the  number  of  terms  with  non-zero  coefficients.  Programs  CADNIP  and 
BADMEP  were  utilized  by  Dorothy  Gillette  (LKB). 

3.3  Results  and  Discussion 

Tables  4 and  5 compare  the  selected  and  the  full  SElll  prediction  errors 
for  the  DB  satellites  and  OVl-17,  respectively.  For  many  of  the  DB  runs  the 
resulting  prediction  error  differences  between  full  and  selected  geopotentials 
lie  within  the  .5  km  target.  The  others  are  all  within  1 km.  The  OVl-17 
prediction  differences  are  considerably  larger  as  are  the  total  errors. 

Hence  in  both  cases  the  prediction  accuracies  of  the  selected  and  full 
geopotentials  are  nearly  identical. 

The  larger  prediction  differences  for  OVl-17  are  due  to  several  possible  effects; 
non-uniform  time  density  of  observations,  variations  in  fit  span,  non-linear 
effects,  and  the  use  of  only  3 orbital  parameters,  combined  with  the  longer 
prediction  time  - fit  span  ratio.  In  particular  data  gaps  of  up  to  3 hours 
occured  in  the  ADCOM  data,  sometimes  at  the  beginning  or  end  of  a fit  span. 

This  would  effectively  shorten  the  fit  span  to  which  the  prediction 
error  effect  is  highly  sensitive  for  larger  prediction  times.  Non- 
linearities  in  the  orbit  relative  to  the  adjustable  parameters  can  also  be 
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Case  # 

Full  SElll 

Selected  SElll 

1 

-3.2 

-2.9 

2 

-3.3 

-2.3 

3 

16.0 

16.1 

4 

0.7 

1.6 

5 

1.3 

1.7 

6 

-3.8 

-3.1 

7 

-0.5 

-0.9 

8 

-1.1 

-0.4 

9 

0.3 

0.3 

10 

-0.2 

0.3 

11 

3.5 

3.7 

12 

-2.5 

-1.5 

Case  # 


Fun  SElll 


Selected  SElll 


1 

30.3 

29.0 

2 

121.8 

112.5 

3 

-5.0 

-2.4 

4 

12.0 

14.5 

5 

-32.4 

-25.9 

6 

14.2 

22.1 

7 

-15.3 

-10.4 

8 

-24.8 

-14.8 

9 

-5.3 

6.7 

10 

-31.6 

-26.6 

11 

2.6 

0.9 

i 

f 

I 


TAble  5.  3 Day  Prediction  Errors  (km) 

Satellite  OVl-17 


expected  to  increase  for  larger  time  periods.  It  was  for  these  reasons 
that  the  relatively  small  (compared  to  overall  prediction  accuracy)  prediction 
difference  criterion  was  used  for  OVl-17.  Alternatively,  one  can  develop  a 
model  to  include  the  above  mentioned  effects.  It  was  not  the  intention  to 
do  so  in  this  work,  but  rather  to  develop  as  simple  a method  as  possible  to 
achieve  a desired  goal. 

The  method  of  geopotential  term  selection  described  here  is  not  put  forth 
as  the  "best"  method  to  be  used  in  all  circumstances.  Certainly  different 
goals  (such  as  accuracy  in  the  orbit  parameters,  or  the  orbit  itself  over  the 
fit  span)  would  suggest  different  criteria.  It  is  believed,  however,  to 
represent  a good  starting  point  and  is  certainly  adequate  for  satellites 
where  errors  from  other  sources  dwarf  those  due  to  the  geopotential.  The 
analysis  presented  in  Section  3.1  should  also  serve  as  a good  starting  point 
for  rapidly  estimating  modeling  error  effects  on  orbit  determination  and 
other  analyses  such  as  atmospheric  density  determination  from  satellite  drag. 
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4.0  Scintillations 
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Problem  No:  4844 

s 

r 

I 

' 4.1  Low  Latitude  Study 

The  system  of  programs  developed  for  creation  and  analyses  of  the  high 
latitude  scintillations  data  base  (References  1 and  2)  were  adapted  and 
^ applied  to  the  equatorial  scintillations  data  base  for  two  stations-- 

; Huancayo  {LES-6  and  ATS-3  satellites)  and  Ghana  (ATS-3  satellite).  The 

( 

! following  programs  comprise  the  scintillations  data  processing  system: 

1.  An  Edited  Raw  Data  file  is  created  using  program  SEQCHK.  The 

raw  data  cards  are  previously  read  in  time  sequence  and  catalogued.  Format 
and  sequence  errors  as  well  as  duplicates  and  invalid  ranges  are  checked. 
Corrections  are  made  over  the  Intercom  System. 

2.  Program  DABASE  creates  the  data  base  with  complete  environmental 


information 

for  each  15-minute  sample. 

Each  sample  was  recorded 

in  BCD 

as  follows: 

IDSTA, 

IDSAT,  IFREQ,  lYR, 

IMO,  IDA, 

UT,  SCIN,  DB,  XKP, 

LF2 

12. 

A2,  13,  12, 

12,  12, 

F5.2,  F4.1,  A2,  F3.1, 

13 

where 

IDSTA: 

Huancayo  = 20, 

Ghana  = 

17 

IDSAT: 

ATS3  = A3, 

LES6  = 

L6 

IFREQ 

= 137  MHz  for  A3, 

254  MHz 

for  L6 

All  OdB  SI 

values  were  discarded.  Ghana 

A3  data  were  converted 

to  dB  values 

by  setting  SI  >.  0.1  and  £ 94.24,  and  then 

SCIN  (dB)  = 15.23077-  /231 .9763-2. 461 54SI 


J 
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Kp  and  2.7  GHz  solar  flux  (LF2)  files  were  provided  by  members  of  the  Space 
Physics  Laboratory. 

Subsequently  all  data  were  compacted  into  one  60-bit  word  with  common 
identifier  characters  entered  in  a leader.  The  format  and  some  details  are 

described  in  Section  4.3,  but  the  programs  below  have  not  yet  been  re- 
vised to  operate  with  data  in  this  format. 

3.  All  further  analyses  of  the  data  were  done  in  local  rather  than 
universal  time.  The  sub-ionospheric  longitudes  used  were 
HA3  = 75®  W,  HL6  = 72.6®  W,  GA3  = 10®  W. 


Program  SIDIST  provides  distributional  breakdowns  in  a form  which  covers 
most  cross-sections  of  interest.  Figure  1 shows  a typical  print-out  for  HA3 
for  month  2,  the  lowest  Kp  range,  and  all  solar  flux  ranges  combined. 
Percentage  occurrences  by  local  time  intervals  and  SI  ranges  are  calculated 
after  sorting  the  data.  The  following  divisions  are  used; 

Months  1,  3(Season),  12  (all) 

Solar  Flux  0 to  95,  96  to  120,  121  up,  all 

Kp  0 to  1-,  1 to  2-,  2 to  3+,  4-  to  5+,  6-  up,  all. 

Local  Time  1 , 3,  24  (al  1 ) 

4.  The  analytical  modeling  effort  requires  that  the  data  be  divided 
into  moderate  size  partitions  in  each  of  which  the  average  or  median  is 
obtained.  Program  SIPART  used  4 Kp  ranges:  0 to  1+,  2-  to  3+,  4-  to  5+,  6-  up 

and  3 S^  ranges:  0 - 95,  96  - 120,  121  up 
for  each  month  and  hour  of  local  time.  Within  each  partition 
S^  and  SI  were  averaged  and  a total  count  was  also  kept  for  possible  future 
weighting.  This  averaged  file  was  saved. 
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Figure  1.  Typical  Print-out  .of  Distributional  Characteristics  c; 
Scintillations  for  Huancayo  ATS3. 
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5.  Program  SIMED  obtained  the  median  scintillation  index  rather  than 
the  average.  This  required  sorting  of  the  data  in  each  partition.  The 
median  and  the  average  were  shown  to  be  similar  for  any  of  the  data  bases, 
and  no  further  studies  were  conducted  using  medians. 

6.  Program  TABSIN  used  the  catalogued  file  to  tabulate  average  or  median 
SI  values  by  partition,  3 months  per  page.  This  concise  format  is  identical 
to  the  output  from  the  modeling  program. 

7.  Program  MODEL  used  the  partition-averaged  file  to  develop  an  ana- 
lytical model  of  scintillation  index  as  a function  of  day  of  the  year,  local 
time,  K and  S..  Out  of  a maximum  of  3456  blocks  (12  months  x 4 K x 3 S. 

p f ' p f 

X 24  local  time)  the  averaged  files  with  at  least  5 observations  in  the 
block  were  as  follows: 

Huancayo  A3  2838  blocks 

Huancayo  L6  2461  blocks 

Ghana  A3  2615  blocks 

The  empty  blocks  generally  occur  in  the  high  Kp  and  occasionally  the  high 
range. 

The  development  and  application  of  program  MODEL  has  been  described 
earlier(‘).  An  additional  feature  now  available  is  a tabulation  of  the 
correlation  coefficients  for  all  the  variables  in  the  search.  Highly 
correlated  variables  can  suggest  reformulation  of  the  model  expression. 

The  final  expression  selected  was  similar  in  form  to  the  high-latitude 
version^^^,  and  Figure  2 summarizes  all  three  cases.  The  standard  deviation 
of  the  residuals  in  each  case  was  as  follows: 

Huancayo  A3  1.77 

Huancayo  L6  1.30 


Ghana  A3 


1.34 


DA  = Day  Number 
HL  = Local  Time  (Hrs . 
As  = Sf/100,  where 


Figure  2.  Equatorial  Scintillation  Index  (dB)  Models 


Since  the  equatorial  scintillations  are  lower  than  those  encountered  in 
the  auroral  regions,  the  lower  residuals  that  were  obtained  were  to  be  ex- 
pected. 

o.  Program  SIPLT  was  used  to  simultaneously  present  the  model  and 
the  averaged  data  which  had  been  used  as  the  basis.  Figures  3-8  show 
typical  plots  for  months  2 and  6 for  the  three  stations. 

4.2  High  Latitude  Data  Extension  and  Analysis 

The  existing  high-latitude  data  base  was  extended  by  including  the  addi- 
tional year  of  data  acquired  for  Narssarssuaq , Goose  Bay  and  Sagamore  Hill 
since  the  original  data  base  was  created.  Programs  SEQCHK  and  DABASH  were 
used;  where  DABASH  was  modified  from  the  original  program  to  include  esti- 
mation of  near-geostationary  satellite  ephemerides  using  inclination 
and  ascending  node  information  from  widely  separated  element  sets.  Sub- 
ionospheric  latitude  and  longitude,  and  geomagnetic  latitude  and  time  , 

are  obtained  using  routine  SILL2  provided  by  the  Space  Physics  Laboratory,  j 

i 

followed  by  C0RRGM2  and  MAGTIM  based  on  Gustafsson's  1970  Revised  Corrected  ! 

Geomagnetic  Coordinate  System.  1 

.1 

Program  SIDIST  was  used  to  provide  distributional  breakdowns  as  were  done  I 

for  the  equatorial  cases  (Figure  1).  | 

The  Narssarssuaq  - Goose  Bay  - Sagamore  Hill  based  global  model  of  Reference  j 

1,  P.  193  which  had  been  derived  from  the  individual  models  for  the  three  j 

j 

stations  was  found  to  possess  a singularity  near  54®  geomagnetic  latitude. 

The  model  was  therefore  corrected  using  a revised  form  which  avoids  this 
problem.  The  revised  model  is  shown  in  Figure  9. 

4.3  Packed  Data  Bases 

The  limited  resolution  of  all  parameter  values  in  the  data  base,  and 
the  extensiveness  of  the  samples  warranted  packing  the  full  data  bases 
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Figure  4.  Average  SI  Data  and  Model  Plots  for  Huancayo  A3  for  June. 


Figure  7.  Average  SI  Data  and  iiodel  Plots  for  Ghana  A3  for  February, 


Figure  8.  Average  SI  Data  and  Model  Plots  for  Ghana  A3  for  June. 


with  high  efficiency.  Each  of  the  bases  for  the  equatorial  study  were 
collected  into  files  with  headers  and  stored  on  one  magnetic  tape.  The 
same  was  done  for  all  the  high  latitude  data. 

All  data  are  in  integer  form,  obtained  by  multiplying  the  orginal  values 
by  the  scale  factors  described  in  the  file  headers  (see  Figures  10  and  11). 
Information  common  to  all  samples  in  the  file,  such  as  station  ID,  sate- 
llite ID,  Frequency,  Preparation  Date,  and  the  number  of  bits  per 
variable  is  also  contained  in  the  headers.  Each  sample  is  described  in 
a CDC  60-bit  word  by  time(UT),  SI(dB),  Kp  and  2.7  GHz  solar  flux.  The 
special  key  is  presently  unassigned  and  contains  zero.  The  high  latitude 
data  requires  an  additional  CDC  word  for  sub-ionospheric  geographic 
latitude  and  longitude,  geomagnetic  latitude  and  geomagnetic  local  time. 

A system  of  routines  PACKUP  (2123/03/02/4937)  for  handling  these  data  is 
described  in  Section  9. 

Periods  covered  and  total  number  of  15-minute  sample  records  that  have 
been  collected  into  the  data  base  are  as  follows: 


DATA  BASES 

SPAN 

#RE CORDS 

Narssarssuaq 

! 

Sep.  68-Jul. 

196388 

Goose  Bay 

Jan.  72-Jul . 76 

120507 

Sagamore  Hill 

Dec.  69-Jan.  76 

165009 

Huancayo  A3 

May  69-Apr.  76 

1644427 

Huancayo  L6 

Dec.  70-Mar.  76 

123715 

Ghana  A3 

Oct.  71-Apr.  76 

154222 

69 


PF  - Geomagnetic  Latitude 

DA  - Day  Number 

HL  - Local  Time  (hrs. ) 

AS  = (2.7  GHz  Solar  Flux)/100. 


SI  = -2.0+1.2CL+6.5(l-.2FD){l+.6(l-.16FD)(l-.45SL)cos(HL  +2  - .4Kp) 
+.05cos(2(HL-.9))+.02cos(3(HL+3.5))}(EXL)  EXA 


u rv«  ,{.2Kp(l+.25CL)(l-.lXD)  + AS(1-.8SL)(1+.45FD)} 
where  EXA  =2 

and  EXL  = 


where 


!CL=cos.((PF-54.4)  ir/25)  CDL  = 1 + .75SL 

SL=sin((PF-54.4)  ir/25)  SDL  = 1-3SL 

PD=cos(DA+16)+.3cos(2(DA-30));  FD=CDLxPD  XD=SDLxPD 


Figure  9.  Revised  Model  of  Scintillation  Based  Upon 
Narssarssuaq,  Goose  Bay,  and  Sagamore  Hill  Data 
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HUANCAYO  L6 

IDSTA  20  IDSAT 

L6  I FREQ  254 

SCINDEX  DB  PREPARED  07/27/77 

FORMAT  PER  RECORD 

VARIABLE 

SCALE 

BITS 

lYR 

1 

9 

I MO 

1 

6 

IDA 

1 

6 

UT 

16 

9 

SI 

10 

9 

KP 

3 

6 

SF2.7 

1 

9 

SP  KEY 

1 

6 

Figure 

10.  Header  for 

Huancayo  L6  Equatorial  Data  Base. 

Tape  EQSIDB  contains  three  files  of  the  equatorial  SI  data  bases 
for  Huancayo  A3,  Huancayo  L6  and  Ghana  A3.  Each  file  starts  with 
the  heater  which  may  be  printed  out  as  80-column  card  images. 
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Figure  11.  Header  for 

Goose  Bay  A3 

Polar  Data  Base 

Tape  PLSIDB  contains  three  files  of  the  polar  SI  data  bases 

for 

Narssarssuaq,  Goose  Bay  and  Sagamore  Hill. 

Each  file  starts 

with 

the  header  which  includes  data 

on  the  station  ID,  satellite 

ID, 

Frequency  and  Preparation  Date, 

as  well  as 

a description  of 

the 

record  format.  These  may  be  printed  out  as 

80-column  card 

images. 
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b.O  Defense  Meteorological  Systems  Project  (Df^SP)  ' 

' Initiator:  R.  C.  Sagalyn 

j 

Project  No:  DMSP  Problem  No:  4836 

Program  SSIE  was  written  to  process  special  ion  and  electron  sensor 
data  from  a series  of  polar  orbiting  satellites  of  the  SAMSO  DMSP  (Defense 
Meteorological  Systems  Project).  The  program  will  be  run  to  process  real  time  j 

data  as  it  is  received  at  the  Air  Weather  Service,  Offutt  Air  Force  Base,  i 

Nebraska,  but  it  has  been  written  and  tested  at  Air  Force  Geophysics  Labora- 
tory. It  has  therefore  been  necessary  to  simulate  the  AWS  Univac  1110  real- 
time system  on  the  AFGL  CDC  6600  computer.  There  were  two  principal  pro- 
gramming tasks;  data  management  and  analysis.  ■ 

5.1  Data  Management 

Preprocessed  and  Processed  Data  Files 
Preprocessed  ion  and  electron  sensor  data  resides  on  a sector  addressable 
circular  file  at  GWC  (addressing  starts  at  sector  0).  The  SSIE  processing 
program  uses  the  preprocessed  file  as  input,  and  constructs  a file  of  similar 
structure  for  output.  The  output  file  contains  the  same  information  as  the  i 

input  file,  but  in  different  order  and  with  the  results  of  the  analysis  ; 

added. 

A sector  is  28  words  long,  and  the  file  size  is  determined  by  the  number 

of  satellite  readouts  (usually  a readout  is  one  orbit's  worth  of  data)  which 

are  to  be  retained  in  the  file  (usually  about  a day's  worth). Both  input  and  ; 

1 

output  files  are  from  20000-25000  sectors  long.  The  following  description  | 

is  for  the  input  file,  but  much  of  it  applies  to  the  sl 

6^  I 


r 

output  file  as  well.  Notes  on  the  output  file  are  included  as  Section 

5.6. 

The  file  is  circular  in  operation  in  that  new  readouts  received  are 
inserted  where  the  last  readout  ended.  Should  a readout  begin  near  the 
end  of  the  file  and  attempt  to  overflow  the  file,  the  data  will  continue 
writing  at  sector  number  4.  (Sectors  0-3  are  reserved  for  the  directory). 

This  scheme  assures  constant  collection  of  the  most  current  data,  since  only 
the  oldest  readouts  are  being  overwritten. 

5.2  General  File  Structure 

Figure  1 shows  the  general  structure  of  the  circular  file.  Figure  1 is  in- 
tended to  depict  the  general  file  structures  and  how  they  might  appear  within 

1 I 

I a circular  file  after  it  has  been  operational  for  some  time.  Figure  2 de-  I 

[ tails  the  contents  of  the  circular  file.  ■ 

I- 

i A.  Directory  - The  main  function  of  the  directory  is  to  maintain  , 

' i 

current  information  concerning  the  status,  position  and  size  of  each  readout  j 

contained  in  the  file.  The  directory  is  located  in  sectors  0-3  (Figure 

2).  The  first  2 words  of  the  112  word  directory  are  used  by  the  processing  ] 

program  to  maintain  file  position  information  only  for  the  most  current  read-  j 

( 

out  received.  The  next  26  words  are  pointers  and  reflect  the  total  file  ; 

^ content.  Each  valid  pointer  word  represents  a readout.  There  are  a larger  | 

number  of  pointers  than  there  are  readouts  that  the  file  could  actually 
contain.  Each  pointer  word  consists  of  3 subdivisions--!,  four  status  bits, 

I 2.  the  sector  number  where  the  readout  begins,  3.  the  number  of  data  sets  i 

^ in  that  readout.  Whenever  a readout  is  overwritten  the  first  status  bit 

for  its  pointer  word  is  set,  and  that  pointer  becomes  in  effect,  invalid. 

The  remainder  of  the  directory  is  unused  at  this  time. 
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Figure  2.  Contents  of  the  Circular  File  (Continued) 
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B.  Readout--A  readout  is  the  variable  quantity  of  data  received 
during  a particular  satellite  data  acquisition  period.  Each  readout  begins 
at  a variable  sector  location  dependent  upon  the  length  of  previous  read- 
outs. The  file  sector  where  a readout  begins  is  given  in  the  directory 
pointer  word.  Each  readout  consists  of  an  information  sector  followed  by  a 
variable  number  of  data  sets.  Readouts  are  placed  in  the  circular  file  in 
the  sequence  in  which  they  are  received  from  the  spacecraft.  A spacecraft 
having  relatively  older  data  on  its  recorder  could  be  readout  after  a 
spacecraft  having  more  recent  data,  thus  the  readouts  are  not  necessarily 
in  time  sequence. 

C.  Information  Block--The  first  four  sectors  of  each  readout  contain 
the  information  block.  The  function  of  this  sector  is  to  maintain  para- 
metric data  peculiar  to  the  readout.  (Figure  2). 

D.  Data  Set--A  data  set  is  defined  in  the  preprocessed  files  to  be 
one  ephemeris  group  plus  60  seconds  of  sensor  data.  Each  data  set  contains 
392  words  and  occupies  14  sectors  (or  420  words  and  15  sectors  for  the 
output  file). 

E.  Ephemeris  Group--The  first  32  words  (60  words  on  output  file)  are 
defined  as  an  ephemeris  group.  The  first  five  words  of  the  ephemeris 
group  contain  position  information  at  the  time  of  the  first  data  message 
following  the  ephemeris  group.  The  second  five  words  contain  position 
information  at  a time  60  seconds  earlier  than  the  time  of  the  first  data 
message  following  the  ephemeris  group.  That  is,  the  information  corresponds 
to  the  time  of  the  data  in  the  last  data  message  in  the  data  set.  This  is 
true  because  the  spacecraft  tape  recorders  are  not  rewound  prior  to  playback. 
The  data  frames,  within  each  readout,  will  be  inserted  into  the  file  in  re- 
verse time  order,  i.e,  the  first  frame  within  each  readout  will  be  the  last 
frame  produced  by  the  sensor.  The  ephemeris  group  for  180  bit/sec  to  216 
bit/sec  sensors  is  contained  in  the  first  28  words  of  the  data  set,  thus  the 


r 
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sensor  data  begins  with  the  first  word  of  the  next  sector  rather  than 
word  5.  The  block  5D  spacecraft  has  precise  attitude  control  which  main- 
tains the  sensor  pointing  accuracy  to  0.01  degree,  thus  pitch,  yaw,  and  roll 
corrections  are  not  necessary. 

F.  Sensor  Data--l-ol lowing  the  ephemeris  group  are  60  one  second 
data  frames.  The  first  word  of  each  data  frame  contains  the  spacecraft 
24  hour  elapsed  time  counter  followed  by  a synch  code  (Figure  2). 

The  least  significant  bit  of  the  time  code  corresponds  to  2 seconds.  The 
synch  code  is  the  following  logical  pattern  =11111000110111010100. 

^ The  time  synch  word  is  followed  by  the  integral  number  of  36-bit  words  ] 

necessary  to  contain  the  total  bit  output  for  the  sensor  when  interrogated  by  | 

! 

the  spacecraft  processor  for  that  second.  The  sensor  data  words  will  not 
be  in  the  same  sequence  in  which  they  were  generated  by  the  sensor.  The  bit 
stream  is  received  in  reverse  order  at  AFGWC,  that  is,  least  significant 
bit  first.  The  data  formatter  (DF)  reverses  the  bits  in  36-bit  groups  (one 
Uflivac  word).  Thus  the  sensor  words  in  the  preprocessed  circular  file  appear 
with  the  most  significant  bit  first  but  in  a different  sequence.  The  data 
word  sequence  for  various  bit  length  words  is  given  in  Figure  3.  Figure  4 
shows  the  contents  of  a 180  bit  data  block. 

5.3  File  Conventions 

A.  If  a data  drop-out,  i.e.  time  discontinuity,  is  encountered  between 
data  frames,  the  remainder,  (if  any)  of  the  60  second  block  of  data  frames 
will  be  zero  filled,  the  ephemeris  groups  associated  with  that  data  set 
will  still  be  60  seconds  apart  as  if  the  data  had  continued  uninterrupted. 

The  next  data  set  generated  will  start  with  the  first  data  frame  following 
the  discontinuity. 
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Figure  3. 


B.  Data  will  be  written  to  the  file  by  the  preprocessor  in  blocks  of 


4 minutes  of  sensor  data  and  associated  ephemeris  groups,  or  56  sectors. 

The  circular  file  will  be  wrapped,  that  is,  cease  writing  at  the  end  of  the 
file  and  begin  overwriting  at  the  beginning  of  the  file,  whenever  there  is 
not  sufficient  room  remaining  in  the  file  for  a complete  four  minute  data 
block  (56  or  64  sectors).  The  remainder  of  the  file,  between  the  last  block 
written  and  the  physical  end  of  the  file,  is  not  initialized  or  set  to 
zeroes.  It  will  contain  whatever  was  left  from  previous  writes. 

C.  If  there  are  insufficient  data  at  the  end  of  a readout  to  com- 
pletely fill  a four  minute  block,  the  remainder  of  the  last  four  minute 
block  will  be  zero  filled. 

5.4  Univac  1110  Conventions 

The  following  information  is  provided  for  those  not  familiar  with  the  Univac 
1110  system. 

Word  length  is  36  bits 

1 sector  is  28  words, mass  storage  is  sector  addressable 
Negative  values  are  ones  complemented 

Mass  storage  I/O  efficiency  is  achieved  in  groups  of  112  words. 

5.5  DMSP  SSIE  Data  Output 

The  SSIE  experiment  provides  a single  digital  NRZ  (nonreturn  to  zero)  data 
signal,  designated  as  (SSIEDAT)  to  the  OLS  (Operational  Linescan  System). 
This  data  signal  is  transferred  in  phase  with  the  OLS  supplied  bit  clock  in 
bursts  of  180  contiguous  bits  at  a bit  rate  of  1000  Hz  + 1 Hz.  Upon  command 
(SSIERED)  from  the  OLS,  the  SSIE  provides  one  180  bit  data  block  per  second 


with  the  least  significant  bit  (LSB)  occurring  first  in  the  first  word. 


The  SSIE  180  bit  data  block  consists  of  twenty,  9 bit  data  words  as  follows: 
7 Samples  Electron  Data  (Words  1,  4,  7,  10,  13,  16,  19) 

7 Samples  Ion  Data  (Words  2,  5,  8,  11,  14,  17,  20) 

6 Samples  Event  Monitor  (Words  3,  6,  9,  12,  15,  18) 

All  data  is  stored  in  SSIE  shift  registers.  Once  per  second  a read  pulse 
(SSIERED)  of  180  ms  duration  from  the  OLS  allows  the  registers  to  shift  out 
data  acquired  during  the  1000  ms  period  prior  to  the  E^  of  the  read  period. 
Word  1 is  shifted  out  first. 
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5.6  Output  File  Notes 

The  processed  file  (output)  from  the  SSIE  analysis  program  fits  the  des- 
cription of  the  DMSP  preprocessed  circular  file,  with  the  following  ex- 
ceptions . 

1.  The  processed  file  is  22016  sectors  long,  addressable  0-22015. 

Each  data  set  is  one  minute  of  data,  but  is  420  words,  or  15  sectors,  long. 
The  one  additional  sector  is  included  in  the  emphemeris  group  (header) 
making  that  group  60  words  in  length,  rather  than  32  as  in  the  preprocessed 
file. 

2.  As  in  the  preprocessed  file,  one  second  of  data  is  contained  in  six 
36  bit  words,  the  first  being  a time-synch  word,  and  the  remaining  five  con- 
taining 20  nine  bit  words  of  data.  In  the  processed  file,  these  words  and 
the  data  sets  data  blocks  are  in  time  order,  but  within  the  blocks  the 

data  are  in  the  same  order  as  in  the  preprocessed  file. 

3.  The  60  word  data  set  header  (ephemeris  group)  for  the  output  file 
is  as  follows; 

First  Sector  of  Minute 


Word 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


Contents 

Latitude,  radians,  0-PI/2,fp. 
Longitude,  " " " 

Altitude,  NM,  tp. 

Julian  Date,  1-366, integer 
Seconds  from  OOOOZ  on  JD, 0-86399 


As  in  words  1-5 


First 
data/time 
in  minute 


Last 

data/time 
in  minute 
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Word  Contents 

Earth  center  coordinates,  time  of  word  5 
T1  X fp 

12  Y fp 

13  Z fp 

Earth  center  coordinates,  time  of  word  10 

14  X fp 

15  Y fp 

16  Z fp 

17  Orbit  anomaly  angle,  r6d,fp.  Time  of  word  5 

II  * II  II  It  ti  II  II  II  2Q 

19  JD  of  latest  sweep  pair  start  time,  integer,  1-366. 

20  Seconds  from  OOOOZ,  integer  0-86399. 
tQ  designate  tii;’.o  at  midpoint  of  the  electron  sv.eep.  If  the  start  cf  an  electron 

sweep  occurs  within  a minute,  the  time  in  words  19-20  will  be  between  the  time  in 
4-5  and  9-10. 

21  Electron  density  at  tg  for  latest  sweep  not  starting  in  current 

_3 

minute,  CM  ,fp. 

22  Electron  density  at  tg  for  sweep  starting  in  current  minute  if  there 

-3 

is  one  - otherwise,  same  value  as  wd  21.  CM  ,fp. 

23  Ratio  of  electron  density  at  tg , to  measured  current  at  sensor  potential 

I 

equal  to  resting  bias  voltage  for  the  same  sweep  as  in  word  21. 

f 

I 24  As  in  word  23  for  same  sweeo  as  in  word  22. 

I 

^ 25  Electron  temperature  at  tg  for  latest  sweep  - current  minute  sweep 

! it  one  starts  in  it,  deg  K,fp. 

26  Vehicle  potential  with  respect  to  plasma  at  tg,  same  sweep  as  in  25, 

volts,  fp. 


Thi?  pdr'riiv.eters  contained  in  the  following  eleven  words  take  their  values  at  to  + 17, 
midway  into  the  ion  sweep  ollowing  the  latest  electron  sweep  - including  the  current 
minute.  There  are  two  coses  a.  2 or  3 species  present  b.  one  specie  present.  In 
case  b there  are  two  ways  to  compute  temperature,  AMU,  and  vehicle  potential  with  re- 
spect to  plasma,  depending  on  how  the  density  is  estimated.  Both  cases  are  carried  - 
with  subscripts  (ss)  1 or  2. 

27  NI  number  cf  species  present,  integer 
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2 or  3 Species 
N(i)  H + density 


Second  Secto''  of  Minute 


1 

2 

3 

4 

5 

6 

7 

8 
9 


N(2)  He^  density 
N(3)  0^  density 

T(])  Average  ion  temp. 
T(2) 

EM(l)  Average  ion  mass 
EM(2) 

Vehicle  pot  wrt  plasma 


One  Specie 
SS  1 density 


SS  2 density 

SS  1 ion  temp. 

SS  2 ion  temp. 

Ion  mass  for  density  N(l) 
Ion  mass  for  density  N(2) 
SS  1 vehicle  pot 
SS  2 vehicle  pot 


Current  at  0 applied  voltage 
The  following,  Scale  height,  is  computed  at  t + 8.5,  midway  between  the  two  sweeps, 

10  Scale  height  for  2/3  specie  case,  or  SS  1 of  1 specie  case  KM. 

11  Scale  height  for  1 specie  case,  SS  2. 

4.  Directory  sector  is  maintained  for  output  file  just  as  for  input  file. 

5.  Information  block  for  a readout  constructed  for  output  file  aS  for  input 
file,  except  that  words  10  and  11  are  interchanged  since  time  order  is  reversed. 

6.  Discontinuities  in  the  data  on  the  output  file  occur  sporadically.  They 
can  be  detected  in  two  ways: 
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d.  tnd  tiii.c  in  header  of  one  minute  does  not  correspond  to  start  time  of  the 
next.  It  may  be  earlier  or  later. 

b.  The  data  words  will  be  all  zeros  at  the  beginning  of  a minute  of  data. 

Conditions  a.  and  b.  may  occur  singly  or  simal taneously . 

Wi'.en  discontinuities  were  detected  on  the  input  preprocessed  file  during  a 
sweep  pair,  the  sweep  pair  was  considered  invalid  and  the  analysis  not  done  for 
that  pair. 

I 

I ■ 

► 

General  Reference:. 


Air  Force  Global  Weather  Control  Data  Acquisition  and  Processino  Special 
Sensor  Data  Processing  Section,  Offut  Air  Force  Base,  Nebraska;  "Snecial 
Sensor  Data  Preprocessed  Circular  File,  Subsystem  Sneci fi cation" , 

July  3,  1976,  Nebraska 
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Section  6,  Plasma  Bulk  Motion 
Author:  S.  T.  Lai 


t-iqure  contents  or  tne  tircuiar  t-iie 


BEST  AVAILABLE  COPY 


6.0  Plasma  Bulk  Motion 


Initiator:  P.  Wildman 

Project  No:  2311  Problem  No:  4701 

Disturbances  and  irregularities  in  the  ionosphere  and  magnetosphere, 
particularly  in  the  polar  and  auroral  regions,  affect  electromagnetic 
wave  communications,  detection,  tracking,  guidance,  and  early-warning 
defense  systems.  Ambient  plasma  bulk  motion  is  an  important  factor 
in  forecasting  ionospheric  and  magnetospheric  plasma  meteorology  and 
in  understanding  the  processes  leading  to  polar  region  irregularities 
(Figure  1).  Such  irregularities  are  already  known  to  cause  the 
scintillations  in  satellite  signals  received  at  ground  stations 
(Reference  1).  This  work  supports  research  projects  at  AFGL  in  situ 
satellite  sensing  of  plasma  bulk  motion  as  well  as  plasma  densities 
and  temperatures.  The  satellites  concerned  are  S3-2  and  S3-3.  Schema- 
tically, the  sensors  are  divided  into  three  groups:  the  planar  ion 
sensor  cluster,  the  non-planar  ion  sensor  cluster,  and  the  spherical 
Langmuir-type  electron  sensor. 

6.1  Data  Processing 

For  data  processing,  three  types  of  data  tapes  are  required,  viz.,  the 
plasma  data  tape,  the  satellite  attitude  (DM)  data  tape,  and  the  satellite 
ephemeris  (ORMAG)  data  tape.  The  structures  of  these  data  are  documented 
in  the  references  (References  2,  3). 
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The  satellite  signals  consist  of  several  patterns  which  have  to  be 
recognized.  They  are  the  TM  applied  voltage  sweeps,  electron  sensor 
real  sweeps,  electron  calibrates,  ion  calibrates,  and  range  switchings. 
The  flow  chart  of  Figure  2 outlines  the  data  processing  and 
computational  schemes. 

6.2  Ionospheric  Plasma  Electron  Properties 
The  technique  used  in  this  satellite  project  for  measurements  of 
ionospheric  plasma  electron  properties  is  to  analyze  the  Langmuir 
probe  responses.  A TM  applied  voltage  sweep  appears  once  every  128 
sec.  in  the  experiment.  This  voltage  is  also  applied  to  every  real 
sweep  of  the  electron  sensor  in  the  following  128  sec.  period.  The 
voltage-current  response  characteristics  of  the  sensor  output  gives 
information  on  the  electron  properties  of  the  plasma  environment. 
Physically,  the  electron  sensor-current  I in  the  Langmuir  probe 
is  proportional  to  the  velocity,  thermal  distribution,  and  density 
of  the  electrons , i .e. 


where  the  square  root  term  is  due  to  the  electron  velocity,  the  expo- 
nential term  is  the  Boltzmann  distribution  function,  N is  electron 
density,  e is  electron  charge,  and  A is  sensor  aperture  area.  There- 
fore, the  slope  of  logg(I)  plotted  versus  voltage  V(t)  should  yield  the 
electron  temperature  T,  and,  once  T is  found,  N can  be  computed. 


Analysis  of  the  real  data  revealed  that  the  sensor-current  I(V)  deviates 
from  Langmuir's  equation  because  of  several  reasons.  They  are:  (1) 
vehicle  potential  is  nonzero  but  an  environmental  dependent  variable; 

(2)  the  effective  aperture  area  is  affected  by  the  plasma  sheath 
surrounding  the  sensor  and  the  vehicle;  (3)  there  often  exist  at  the 
same  region  more  than  one  species  of  electrons,  characterized  by  their 
densities  and  temperatures;  (4)  geometry  of  grids  in  the  sensor; 

(5)  secondary  electron  emissions  and  also  photo  electron  emissions 
from  the  sensor. 

The  vehicle  potential  variations  in  the  ionosphere  are  related  to  the 
phenomenon  of  spacecraft  charging  (Reference  4).  The  plasma  sheath  is 
a result  of  the  well  known  Debye  shielding  property  of  plasmas. 

The  presence  of  multi -electron  species  is  due  to  electron  flows  in  the 
ionosphere,  especially  in  events  such  as  energetic  electron  precipita- 
tions in  auroral  and  red  arc  regions  and  in  polar  plasma  flow  regions. 
Photoemissions  and  secondary  emissions  are  due  to  solar  photons  and 
energetic  electrons  respectively. 

Mathematically,  the  sensor  current  becomes  expressed  as  the  sum  of 
contributions  of  I-j  and  I2  from  the  two  species  and  from  the  secondary 
and  photoemissions. 

I(V)  = I^(N^,T^,V)  + I2(N2,T2,V)  +I3(V) 


"1 


where  V = + <J>,  where  $ is  the  vehicle  potential,  and  is  the  applied 

voltage  on  the  sensor,  r is  sheath  radius,  R grid  radius,  5 transmission 
factor,  k Boltzmann  constant,  e electron  charge,  N electron  density,  T 
electron  density,  and  i labels  the  species.  Each  sweep  pattern  is  the 
superposition  of  the  sweep  patterns  corresponding  to  two  single  species. 

The  expression  I(V)  is  invariant  under  permutation  of  the  indices  i (=1,2). 
To  determine  the  parameters  that  characterize  the  physical  contents  of 
the  sweeps,  the  method  of  steepest  descent  of  Fletcher  and  Powell 
(Reference  5)  is  used  for  fitting  the  sweep  patterns.  The  parameters  to 
be  determined  are  , T^. , <!>  , r,  a and  b.  Occasionally,  the  results 
of  a sweep  have  to  be  rejected  because  of  over  saturation  of  data, 
broken  sweep  due  to  existence  of  time  gao  in  which  no  signal  was  received 
on  the  ground,  too  much  noise  in  data,  or  simply  overshooting  in  the 
rapid  descent  technique.  Usually,  the  sweep  results  are  very  useful. 

The  scheme  of  electron  density  correction  computation  scheme  at  the  end 
of  every  real  sweep  is  shown  in  Figure  3. 

6.3  Determination  of  the  Plasma  Flow  Direction 
In  this  section  we  derive  expressions  for  the  pitch  andyaw  of  the 
plasma  flow  relative  to  the  spacecraft  coordinate  system,  using  the 
currents  flowing  to  an  array  of  the  positive  ion  sensors.  The  space- 
craft velocity  relative  to  the  plasma  contains  components  due  to  the 
spacecraft  velocity  and  to  any  drift  or  'bulk'  motion  of  the  olasma, 
(References  8,  9).  To  evaluate  the  component  due  to  plasma  motion  it  is 

necessary  to  know  the  spacecraft  attitude  and  velocity  in  an  external 
coordinate  system,  and  this  is  dealt  with  in  Section  6.4.  First,  however, 
we  must  derive  the  pitch  andyaw  angles  of  the  plasma  flow  in  the  space- 
craft geometric  coordinates. 
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Electron  Sensor  Data 
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In  Figure  4,  we  show  x y z,  the  spacecraft  coordinates;  ON  is  the  normal 
to  the  planar  sensor  aperture,  and  OQ  (q)  is  the  relative  plasma  flow 
vector  with  pitch  (E^  + P^)  and  Yaw  Yq'  Angle  is  a time-dependent 
measure  of  the  rotational  position  of  the  spacecraft  as  it  spins  about 
the  y axis.  OV  (v)  is  the  satellite  velocity  vector,  which  at  time 
t is  inclined  to  the  x z plane  at  angle 

For  unambiguous  measurement  of  Pq  and  two  orthogonal  pairs  of  such 
I sensors  are  required  and  the  actual  configuration  used  is  shown  in 

Figure  5. 

The  lower  three  sensors,  #1,  #2,  and  #3  measure  yaw,  and  the  vertical 
pair,  #2  and  #4  measure  pitch  as  the  array  faces  the  plasma  flow  vector 
once  during  each  spacecraft  spin. 

If  the  centerline  of  the  array  (line  bisecting  normals  of  sensors  #2 
and  #4),  nominally  on  the  x axis,  is  rotated  y degrees  about  the  y axis 
and  c degrees  about  the  z axis,  then  we  have  fully  described  the  alignment 
of  the  sensor  array  relative  to  the  spacecraft  axes.  The  angle  between 
the  normals  of  any  two  adjacent  sensor  is  2a  so  that  the  values  of 
P^  and  Y^  for  the  four  sensors  are  as  follows: 

Sensor  No.  1 

Pg  -(a-y) 

Yj  -(2a+e) 

TABLE  1 . Skew  Sensor  Pitch  and  Yaw  Angles 


2 3 4 

-(ot-y)  -(a-y)  (a+y) 

-€  -(2a-e)  -e 
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Considering  again  the  ion  flow  to  a single  sensor,  we  assume  that  the 
mean  thermal  (isotropic)  velocity  of  the  ions  is  much  less  than  the 
spacecraft  velocity  (Sagalyn  and  Smiddy,  Reference  7)  so  that 
I,  the  sensor  current  is  given  by: 

I = C Ne  Aq  Cos  <1) 

where: 

C = transmission  of  the  aperture 

N = ambient  ion  density 

e = electron  charge 

q = relative  velocity  between  plasma  and  spacecraft 
A = aperture  area 

p = angle  between  q and  aperture  normal  (NOQ  in  Figure  4) 


Using  the  direction  Cosines  of  OQ  and  ON  to  derive  Cos  <J),  we  obtain: 

I = 5 Ne  Aq  [cos  Cos  Cos(Pg  - P5  + E^)  + Sin  Y^  Sin  Y^] 
Let  n = Ne  Aq 

X » Pji  t a - Y + Ej 

and  Pg,  Yg  are  the  Pitch  and  Yaw  of  the  sensor  normal.  From  the  angles 
given  in  Table  l we  obtain  the  following  relations: 

■ (I,  - * Ij) 

l-j  1 2 ~ ^ ^ ^ Cos  2 ot 


Now  let  tan  <I>  = tan  Y/Cos  X 
then: 

R21  = tan  a tan  (a  + e + 4') 
R^^  = tan  a tan  (a  - £ - <I>) 
R-j^  = -tan  2 a tan  (e  + <I>) 
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so  that 

R,/,  r-T— B — = tan  a tan  (a  - X) 
I + X *'24 


-R94  1 - X ' 

tan  (Pq  + - y)  = V + X R24  I 

The  choice  between  R21  and  R^^  depends  on  which  side  of  the  center- 
line  the  plasma  flow  is  located. 

6.4  Derivation  of  Plasma  Velocity  Vector 

In  Figure  6 we  show  the  vector  triangle  containing  the  true  plasma  flow 
vector  f together  with  v and  q from  Section  6.2.  The  angle  6,  between 
V and  q is  determined  knowing  Pg  and  Yg  from  Section  6.2. 

In  Figure  6 we  now  have: 

f2  = q2  + yZ  _ 2q  V (Cos  Cos  Yg  Cos  Pg  + Sin  Sin  Yg) 

The  direction  cosines  of  the  true  plasma  flow  are  given  by: 

Cos  a = |y  Cos  Cos  - q Cos  Yg  Cos  (E^  + Pg)j  /f 

cos  e = ^V  sin  - q sin  Ygj  /f 

cos  u = jv  cos  Sin  - q cos  Yg  sin  (E^  + Pg)j  /f 
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The  pitch  Pp  and  yaw  Yp  angles  of  the  true  plasma  flow  are  related  to  these 
direction  cosines  by: 


cos  Yp  cos  (E^  + Pp)  = cos  a 

sin  Pp  = cos  6 

cos  Yp  sin  (E^  + Pp)  = cos  y 

Thus,  the  pitch  and  yaw  angles  can  be  determined  in  terms  of  the 
magnitudes  and  directions  of  relative  plasma  flow  and  vehicle  velocity; 


tan  (E^  + Pp)  = 


V cos  0^  sin  - q cos 

V cos  0^  cos  - q cos 


V„  sir  (£j  + P ) 


Yq  cos 


tan  Yp  = 


V sin  - q sin  Yj. 


cos  + q cos  Yg  - 2Vq  cos  cos  Pg  cos  Y^ 


6.5  Blockage  of  Ion  Flows  and  Non-linear  Function  Fitting 
After  the  careful  analysis  of  a considerable  number  of  orbits  of  real 
satellite  data,  it  became  clear  that  for  cases  where  the  Mach  number 
of  the  satellite's  motion  relative  to  the  ambient  plasma  is  large 
(M  »1),  the  current  measured  by  a sensor  is  severely  limited  when  the 
plasma  flow  vector  intercepts  a portion  of  the  spacecraft  surface 
reaching  the  surface  aperture.  When  the  flow  vector  does  not  intercept 
the  satellite  surface,  the  current  does  behave  similarly  to  that  detected 
by  an  "unobstructed"  sensor,  as  formulated  by  Lai  et  al  (Reference  9). 
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For  the  case  when  the  Mach  number  of  the  satellite  motion  relative  to  the 
ambient  plasma  is  near  or  greater  than  1 (M  ^ 1),  the  current  measured  by  a 

sensor  is  less  than  that  which  would  be  measured  by  an  "unobstructed"  sensor 
for  all  directions  of  the  flow  vector. 

To  model  exactly  the  shapes  of  the  obstructions  of  flow  to  the  sensor 
would  require  an  algorithm  that  would  not  be  efficient  for  the  analysis 
of  a large  volume  of  data.  Instead,  the  obstructions  are  approximated 
by  a series  of  wedges  centered  on  the  sensor's  apertures  (Reference 
10).  The  arrangements  of  the  model  wedges  obstructing  the  fields  of 
view  of  sensors  are  shown  in  Figures  7 and  8.  As  a result,  the 
analytical  expression  of  the  current  measured  by  a sensor  looking  out 
of  the  satellite  spin  plane  becomes  (Reference  10); 


lout  ” NeA(a/ff)5y v {trfU'a  (v.^  cot  ± q cos  rj)) 


+ erf  [vTa  (v^  cot  <!>2  - ^ ?)]  } exp  - q cos  $)]  • 


Similarly,  the  analytical  expression  for  the  in-plane  sensor  current 
becomes 

1 

Ijp  “ NeA  J"  {crf[\/a  cot  + q cos  n)) 

o 

+ erf  (>/a  (i/^  cot  - q cos  5))  ) exp  (-a  (i/  - q cos  e)^]  dv 

z z 
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These  integrals  can  be  calculated  with  good  accuracy  with  Gaussian 
Quadrature.  The  technique  to  determine  the  plasma  flow  parameters  is  to 
use  non-linear  function  data  curve  fitting  methods.  The  parameters 
are  found  by  a search  in  the  vector  space  represented  by  them  for  the 
"solution"  that  yields  a minimum  in  the  difference  between  the  expected 
and  real  values  of  the  function.  The  search  algorithm  is  based  on  the 
method  developed  by  Fletcher  and  Powell  (References  5,  6).  Using  the 
methods  described  in  this  section,  we  have  computed  the  plasma  flow 
data  in  many  orbits  (Figures  9,  10).  In  particular,  a very  high  (9.8  km/s.) 
flow  velocity  of  plasma  was  discovered  to  be  associated  with  a large 
poleward  electric  field  structure  (see  Section  7)  at  the  plasmapause 
during  the  May  1976  magnetic  storm  (Reference  11). 
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Figure  10.  Plot  of  Geophysical  Parameters  Computed  from  Electron  Sensor  Data 
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7.0  Ionospheric  Electric  Fields 

Initiator:  M.  Smiddy 

Project  No:  2311  Problem  No:  4702 


The  reconnection  of  the  magnetic  field  lines  frozen  into  the  solar  wind  with 
the  earth's  field  lines  at  the  magnetopause  leads  to  Urge  scale  induced 
electric  fields  which  are  subsequently  mapped  down  into  the  polar  ionosphere. 
In  order  to  measure  the  global  distribution  of  these  fields  and  their 
associated  currents,  the  Air  Force  Satellite  S3-2  was  launched  into  a nearly 
polar  orbit  late  in  1975. 

7.1  Rotating  Satellite  Boom  Sensor  System 

The  satellite  featured  three  pairs  of  deployable  booms,  each  pair  being 
orthogonal  to  the  others.  Carbon  coated  sensor  spheres  were  mounted  on 
the  ends  of  the  booms.  Data  were  obtained  from  the  axial  pair  and  from  one 
of  the  wire  boom  pairs  in  the  satellite  spin  plane.  The  frequency  of  the 
data  from  the  sensor  was  predominantly  spin  frequency.  However,  other 
frequencies  corresponding  to  various  oscillation  modes  may  be  present, 
especially  during  the  transient  response  periods  following  boom  deployment, 
spin  orientation  maneuver,  etc.  A comprehensive  study  of  the  mode  dynamics 
of  the  rotating  satellite  was  carried  out  in  order  to  analyze  all  mode 
frequencies  of  the  system. 

Lagrangian  equations  of  motion  and  Laplace  transform  techniques  are  used 
for  computations  of  the  transient  response  of  the  satellite  system. 
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For  full  digital  simulation  of  satellite  dynamics  including  nonlinear  | 

mode-mode  interactions  and  combinations  of  various  boom  deployment,  a set  | 

of  Lagrangian  equations  is  solved:  ! 


cL 

dt 


L ({x^},  a,.}  ) = F ( {x. },{;.}) 


where  i,j  = 1,  . . . N.  and  {x^-}  is  a set  of  N dynamical  variables  for  a 
spinning  satellite  system. 


Numerical  solutions  of  these  coupled  nonlinear  differential  equations 
by  means  of  Hamming's  modified  predictor-corrector  method  (Reference  1) 
provide  the  time-dependent  development  of  the  mode  dynamics  of  the  rotating 
satellite  under  specified  boundary  conditions.  This  technique  is  useful 
for  detecting  any  mode  frequency,  instability  side  drift,  excessive  precession, 
and  damping  characteristics  of  satellite-boom  systems  (References  2,3,4). 

7.2  Effect  of  Earth's  Rotation 

The  electric  field  induced  in  the  dipole  sensors  of  a satellite  is  of  sine 
wave  form  because  of  the  spin  of  the  satellite.  In  the  regions  of  low 
latitude  (equatorial)  and  low  altitude,  it  is  well  known  that  the  electric 
field  is  small  with  lack  of  significant  plasma  motion,  so  that  the  predominant 
electric  field  is  due  to  the  Earth's  rotation.  Two  important  remarks  should 
be  mentioned: 

1.  The  electric  field  £ measured  in  a moving  frame  is  related 
to  the  electric  field  Eq  in  a rest  frame  by  (Reference  5). 
i = Eo  + V X B 


j 
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where  ^ is  the  magnetic  field  in  the  rest  frame,  and  is  the  velocity  which  is 

small  (V  « C)  otherwise  there  is  an  additional  relativistic  factor  in  the 
denominator. 


r 


t 

j 2.  The  velocity  Vp  in  a rotating  frame  r is  related  to  the 

velocity  in  a fixed  frame  by  (Reference  6). 

Vp  = V - X r 

where  ^ is  the  rotation  angular  velocity  and  £ is  the  position  vector. 

Thus,  the  induced  electric  field  due  to  the  Earth's  rotation  is 
E_  = ^ X B 

This  induced  field  £ is  always  subtracted  out  in  order  to  bring  cut  the 
genuine  electric  field  from  the  processed  result. 

7.3  Program  System  for  Data  Reduction  and  Analysis 

For  data  reduction  and  computation,  three  types  of  input  data  are  required. 
They  are  1)  electric  field  experiment  data,  2)  satellite  sensor  attitude 
(OM)  data,  and  3)  ephemeri s/magnetic  field  (ORMAG)  data.  The  structures 
of  these  data  are  tabulated  in  the  references  (References  7,8,9). 

The  electric  field  experiment  data  are  in  the  fifth  file  of  every  tape  of  the 
satellite's  raw  data.  The  experiment  data  contain  six  channels  of  sensor 
signals  with  various  high  and  low  gains.  These  signals  are  given  in  telemetry 
counts,  which  can  be  converted  to  voltages  by  a simple  conversion  factor. 

One  unusual  feature  in  the  S3-2  electric  field  data  is  the  appearance  of 
ion  (especially  oxygen  ion)  induced  DC-offsets  detected  by  the  axial  booms. 
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Such  a feature  has  never  been  detected  previously,  because  of  the  relatively  low 
sensitivity  of  the  instruments  on  previous  satellites.  The  DC- 


offsets  vary  with  time  and  are  due  to  variation  of  the  contact  potentials 
on  the  sensors  in  the  ionized  environment  (Figure  1).  No  constant  direct 
correlation  exists  between  the  offsets  and  spin  plane  dipole  measurements 
(Figure  2). 

For  lower  latitudes  (<  60°),  where  the  earth's  magnetic  field  lines  are 
closed,  the  ambient  electric  field  is  vanishingly  small.  Therefore,  any 
voltage  that  remains  after  subtraction  of  the  induction  field  (^  x can 
be  attributed  to  the  spheres'  contact  potential.  By  deleting  regions 
above  and  below  60°,  and  by  fitting  a slowly  varying  function  to  this 
voltage,  one  can  determine,  approximately,  the  contact  potential  for  an 
entire  orbit. 

The  entire  analysis  is  presently  executed  by  means  of  a highly  integrated 
system  of  programs  as  depicted  in  the  flow  chart  of  Figure  3.  The  functional 
form  is  generated  for  the  contact  potential.  This  function  is  then  used 
in  the  MAIN  processing  system. 

r 7.4  Program  Description 

Initial  processing  converts  the  measured  voltage  difference  between  spheres 
, by  the  appropriate  amplifier  gain,  removes  amplifier  biases,  subtracts  the 

induction  field,  and  eliminates  calibration  sequences. 
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Figure  1.  Channel  6 Data  of  The  noisy  regions  are  due  to  Auroral  Electric  Fields. 
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Figure  3.  Flow  Chart  of  Program  System  for  Electric  Field  Analysis 
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Interactive  processing  displays  the  intermediate  data  on  a CRT  graphics 
console,  permits  editing,  and  performs  various  functional  fits  to  the 
edited  data  while  allowing  the  user  to  vary  pertinent  parameters.  Plots 
of  the  final  fit  and  a deck  of  coefficient  cards  are  output  on  request. 

The  program  also  has  several  filtering  routines  which  can  be  used  for 
correlation  analyses.  A non-interactive  version  can  also  be  used  optionally. 

Since  one  set  of  sensors  in  the  spin  plane  is  inoperative,  an  instantaneous 
determination  of  £(t)  is  not  possible.  However,  if  £(t)  varies  at  a rate 
slower  than  the  spin  frequency  of  the  satellite,  then  the  field  can  be  i 

resolved.  More  specifically,  if  the  frequency  spectrum  of  E(t),  defined  by 

is  bounded  so  thatll(w.)  1 = 0 for  oo  > w , then  one  can  determine  E 
— s — 

unambiguously.  For  most  orbits  studied,  this  seems  to  be  the  case.  The 
Main  Program  System  then  determines  £ by  fitting  a function  of  the  form 
V(t)  = A(t)  sin  (ojgt)  + B(t)  costjgt)  C(t) 
to  the  spinning  boom  voltage.  A(t)and  B(t)  are  expressed  1n  terms  of  a 
predetermined  number  of  Chebyshev  polynomials,  such  that  their  frequency 
content  does  not  exceed  Since  the  voltage  measured  is  given  by 

V(t)  = £(t)  ‘^(t)  , where 

R(t)  = Rq  (cos  (Ugt)  R + sin(a)gt)  y), 

one  sees  that  A(t)  and  B(t)  are  the  x and  y components  of  £,  expressed 
in  a coordinate  system  defined  at  t=o.  The  third  component  of  the  electric 
field,  is  measured  directly  by  the  axial  boom. 
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After  subtraction  of  (he  induction  field  (V  x B)  and  determination  of  the 
empirical  offset,  resultant  vectors  are  transformed  to  various  physical 
coordinate  systems,  written  on  a disc  file,  and  plotted  as  strip 
charts  and  polar  plots.  In  the  strip  charts,  the  ionospheric 
electric  field  magnitude  and  components  are  plotted  for  a whole  orbit  in 
Earth  centered  inertial  geomagnetic,  satellite  local  vertical,  or  optionally 
solar  magnetic  coordinates.  A typical  plot  is  shown  in  Figure  4.  In  the 
polar  plots,  the  electric  field  geomagnetic  components  in  the  polar  caps 
are  mapped  on  a background  of  invariant  latitude  and  magnetic  local  time 
(Figure  5,  6).  For  dawn-dusk  orbits,  the  integration /E.^t  of  the  polar 
cap  electric  fields  is  also  plotted  (Figure  7). 

7.5  Results  and  Discovery 

Using  the  techniques  outlined  in  this  section,  many  orbits  of  satellite 
data  have  been  processed  and  the  electric  fields  mapped.  In  particular, 
the  largest  electric  field  structure  ever  detected  at  low  altitudes  in  the 
ionosphere  has  been  discovered  (References  10,  11).  The  event  took  place 
during  the  May  1976  magnetic  storm  and  the  field  structure  was  in  the  evening 
sector  near  the  ionospheric  projection  of  the  plasmapause  (Figure  8,  9). 

The  largest  pulse  was  280  mV/m,  directed  toward  magnetic  north  at  an  altitude 
of  1460  km  with  a hundred  kilometer  latitudinal  extent.  A very  high 
drift  velocity  (about  9.8  km/sec)  of  the  plasma  flow  (see  Sec.  7)  was 
associated  with  the  pulse,  and  the  flow  directed  westward  along  L shells. 
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Figure  8.  Electric  Field  Pulse  during  the  May  1976  Magnetic  Storm.  The  smoothly  varying  sine 
wave  is  primarily  due  to  the  field. 
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8.0  Mul ti -Spectral  Measurements  Program 


Initiator;  E.  Robinson 

Project  No;  2123  Problem  Mo:  4937 

This  program  consists  of  a series  of  rocket  flights  in  which  the  payload 

separates  into  2 stages  - a sensor  module  which  will  contain  sensing 

instrumentation  and  a target  module  which  the  sensing  instruments  will 

be  viewing.  There  will  be  six  radars  tracking  the  vehicle  - 3 for  i 

i 

each  module.  An  attitude  control  system  will  orient  the  sensor  to  ' 

j 

point  at  the  target.  Each  flight  will  last  approximately  600  seconds.  | 

8.1  System  Studies 

The  data  reduction  system  designed  according  to  SUA  specifications  is 

shown  schematically  in  Figure  1.  After  decormutati on  the  input  is  in  the 

form  of  digital  tapes,  2 from  the  sensor  module  (TLl,  TL2)  and  1 from 

the  target  module.  i 

i 

The  data  must  be  unpacked  and  separated  to  form  a number  of  individual 
raw  data  bases  for  further  processing.  Individual  raw  data  bases  are 
then  further  reduced  through  demodulation  and  filtering,  background  and/or 
dark  count  determination  and  subtraction,  and  application  of  pre- 
flight calibration  data.  Attitude,  trajectory,  and  target  data  are  merged 
with  sensor  data  in  a single  file  for  correlative  analysis  (i.e.,  raw  i 

or  reduced  sensor  data  vs.  altitude).  Geophysical  environmental  data,  | 
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such  ds  solar,  lunar  and  stellar  positions  may  also  be  used  in  the  data 
analysis,  although  not  necessarily  merged  into  the  final  data  base. 
Section  8.3  describes  the  preliminary  file  design. 

In  the  initial  phases  of  the  data  reduction  a quick-look  study  is 
carried  out  covering  the  raw  data  before  calibrations  and  using 
preliminary  attitude  and  trajectory  information.  Sample  plots  and 
printouts  are  generated  for  review  by  the  experimenters.  This  phase 
is  expected  to  last  20  days. 

8.1.1  Ultraviolet  Sensors 

The  UV  instrumentation  will  consist  of  3 photometers  to  nrovide  spatial 
resolution  and  a spectrometer  for  spectral  resolution. 

8. 1.1.1  Photometers 

There  will  be  two  medium-sensitivity  photometers  (MSP)  and  one  high- 
sensitivity  photometer  (HSP).  Each  will  contain  a square  array  (10x10) 
of  100  picture  elements  (pixels)  each  with  a field-of-view  of  .2°x.2° 
and  a filter  wheel  of  4 filters.  Readings  will  be  .1  sec  sample-and- 
hold. 

8. 1.1. 2 Spectrometer 

The  spectrometer  consists  of  a linear  array  of  36  spectral  elements 
(specels)  each  responding  to  a different  wavelength  band  depending 
on  the  position  of  the  source  in  the  1.5°xl.5°  field  of  view.  As  with 
the  photometers  the  data  is  .1  sec  sample-and-hold. 


8.1.2  IR  Sensors 


The  first  MSMP  flight  payload  will  contain  4 infrared  instruments:  a 
full-field  radiometer  to  measure  radiation  in  2 wavelength  regions 
(SWIR  and  MWIR),  2 CVF  spectrometers  to  measure  distribution  from  2 to 
5 pm  and  from  7-22  ym,  and  a spatial  radiometer  to  provide  a soatial 
map  of  both  SWIR  and  MWIR. 

8. 1.2.1  Full-Field  Radiometer 

The  input  from  the  full-field  radiometer  consists  of  4 prime  data 
channels  (a  high  gain  and  a low  gain  signal  for  each  of  the  two  wave- 
length regions)  plus  a 50  Hz  reference  sine  wave  signal,  each  sampled 
at  1,000  Hz.  The  data  signals  consist  of  ampli tude-modulated  signals 
at  a carrier  frequency  of  50  Hz.  These  are  to  be  input  into  the 
demodulation/filtering  software  developed  by  the  Optical  Physics 
Division  of  AFGL  and  shown  schematically  in  Figure  2^^^.  This  software 
package  has  been  incorporated  into  the  full  data  reduction  system. 
Briefly,  the  output  of  the  multipliers  contains  a d.c.  component, 
which  contains  the  primary  useful  information,  plus  comoonents  at 
multiples  of  the  carrier  frequency  (50  Hz).  The  higher  frequencies 
are  subsequently  filtered  out  with  the  low-pass  filters. 

8. 1.2. 2 CVF  Spectrometers 

The  input  of  each  spectrometer  consists  of  high-gain  and  low-gain  data 
channels  plus  an  amplitude-modulated  square  wave  spectral  position 
reference  signal  (Fig.  3^^^).  The  fall-off  from  a 10  volt  peak  in  both 
data  and  reference  channel  marks  the  beginning  of  each  new  spectral 
scan.  Pre-flight  calibrations  will  provide  information  on  the  wavelength 
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Figure  2.  Demodulation  of  Radiometer  Signal 


Figure  3.  CVF  Spectrometer  Data  Channels 


for  each  positive  and  neqative  transition  in  a scan.  All  signals  are 
sampled  at  2,000  Hz. 


8. 1.2. 3 Spatial  Radiometer 

The  spatial  radiometer  measures  the  spatial  distribution  in  the  SWIR  and 
MWIR  regions.  This  is  done  with  4 vertical  arrays  of  10  detectors  each, 
2 arrays  per  wavelength  region,  designated  as  small  (.2°x.2°  field-of- 
view)  and  large  (.4°x.4°).  The  detectors  sweep  horizontally  over  a 
6°  range,  1 sweep/sec,  250  data  points/sweep  for  each  detector.  This  is 
the  only  instrument  whose  data  is  to  be  transmitted  via  TL2.  All  other 
sensor  data  is  via  TLl . 

8.1.3  Attitude/Trajectory  Data 

The  following  parameters  are  required  to  define  the  attitude  and 
trajectory  for  full  analysis  of  the  data.  The  appropriate  attitude/tra- 
jectory raw  data  are  to  be  processed  separately  to  produce  a data  base 
with  the  required  parameters  to  be  merged  with  the  sensor.  The 
trajectory  will  be  processed  by  the  AFGL  trajectory  determination 
system,  possibly  including  programs  DRIVEA,  DRIVEB,  and  DRIVEC.^^^ 

The  appropriate  raw  attitude  parameters  are  to  be  taken  from  the  sensor 
and  target  telemetry  tapes  and  processed  by  SUA's  attitude  determination 
group  to  produce  a tape  containing  the  required  attitude  parameters. 

The  trajectories  are  also  required  for  this  phase. 

The  principal  feature  of  the  attitude  requirements  is  the  sensor 
coordinate  system  (Figs.  4 and  5).  The  z'  axis  corresponds  to  the 
line-of-sight  for  the  instrument  while  the  x'  and  y'  axes  define  the 
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Figure  4.  Transformation  from  ECI(xyz) 
coordinates  to  sensor  (x'y'z') 
coordinates . 


Figure  5.  Thrust  Vector  in  Sensor 
Coordinate  System. 


orientation  of  the  sensor  plane.  It  will  be  necessary  to  define  a 
coordinate  system  for  each  instrument.  This  will  require  alignment 
data  relative  to  the  sensor  module's  body  axis,  including  definition 
of  the  x'  and  y'  axes.  Knowledge  of  the  sensor  coordinate  system  permits 
one  to  know  where  the  target  and  stars  are  in  the  sensor's  field  of 
view  and  thus  facilitates  interpretation  of  the  data. 

AHITUDE  DATA 

(1)  Angles  relating  sensor  coordinate  system  x'y'z'  (z'  = line  of 

sight)  to  ECI  system  xyz  (Figure  4^^^)  (0.2°  accuracy) 

i 

1. 

I a)  Right  ascension  of  line  of  sight  = 4)  - ^ 

b)  Declination  of  line  of  sight  = ^ - 0 

c)  Rotation  angle  of  x'  axis  from  line  of  nodes  = \l) 


(2)  Thrust  vector  direction  (a, 8)  in  sensor  coordinate  system  (Figure  5) 
and  roll  angle  of  target  (1°  accuracy) 

(3)  Target  nozzle  position  (A.y)  (Figure  6)  in  sensor  plane  (0.2° 
accuracy) 

(4)  Earth  limb  view  angles  (C.2°  accuracy) 

a)  Angles  relating  sensor  coordinate  system  to  local  vertical  (z)  - 
east  (x)  - north  (y)  coordinate  system  (Figure  4) 

1)  View  azimuth  from  true  north  = it -4) 
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1 


R X TARGET 

I 

I 

I R y 


Figure  6.  Target  Position  in  Sensor  Plane. 
R=Target-Sensor  Slant  Range. 
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1 

2)  Elevation  rel.  to  horizontal  = ^ - 6 

3)  Skew  angle  = ip 

I 

b)  Angle  (6)  between  los  and  hard  earth  tangent  (Figure  7) 

c)  Tangent  height  (BT')  (Figure  7) 

d)  Distance  (OB)  along  los  from  sensor  to  point  nearest  earth 
surface  (Figure  7) 

(5)  Angle  between  target  thrust  and  target  velocity. 


0 HORIZONTAL 


Figure  7.  Earth's  Lipib  View  Angles 
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TRAJECTORY  DATA 


(1)  Position  and  velocity  of  both  vehicles  in  ECI  coordinates 


(2)  Target-sensor  slant  angle 


(3)  Altitude  of  both  vehicles 
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8.2  Calibration  Data  Processing 

Data  processing  support  is  necessary  for  some  pre-flight  calibrations. 
These  calibrations  can  be  conveniently  conducted  using  the  existing 
telemetry  system  along  with  a few  manually  entered  values  indicating 
pertinent  experiment  parameters.  Programs  to  summarize  and  present 
the  results  of  calibration  tests  are  required,  and  to  some  extent,  are 
similar  to  the  programs  used  in  the  flight  data  reduction. 

8.2.1  UV  Photometers 

Pre-flight  calibrations  have  been  carried  out  at  AFGL  for  the  medium 
sensitivity  photometers  MSP-1  and  MSP-3.  Figures  8 and  9^^^  show  the 
data  processing  plan.  In  each  experiment  a beam  of  known  wavelength 
is  mounted  on  a rotating  gimbal  so  that  it  scans  the  entire  array  of 
pixels.  During  a scan  the  gimbal  stops  when  the  beam  points  at  the 
center  of  a pixel,  defining  the  stationary  gimbal  data  to  be  used  to 
determine  its  peak  sensitivity.  Cross-talk  (response  of  pixels  not 
illuminated)  can  also  be  noted.  Moving  gimbal  data  permits  measure- 
ment of  each  pixel's  field  of  view. 

8. 2. 1.1  Program  PHCLFG 

PHCLFG  (Photometer  Calibration,  Fixed  Gimbal)  was  written  to  provide 
calibration  data  for  the  UV  photometers  associated  with  MSMP. 
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N^j(a,AX,I) 


RAW 

FILE 


STATIONARY 

GIMBAL 


MOVING 

GIMBAL 


LISTINGS : 

STANDARD  DEV. 


DATA  BASE : > 


PLOTS,  LISTINGS,  DATA  BASE: 


Sij  (a,AX,I) 


Figure  8.  General  flow  of  UV  Photometer  Calibrations 


Functional  Description 

Dark  Count:  Initially,  many  frames  of  dark  count  are  taken.  The  avei^age 
and  standard  deviation  for  each  pixel  is  calculated  and  printed  out. 

Normal  Data:  At  each  fixed  gimbal  position,  the  average  count  (dark 
count  subtracted)  and  standard  deviation  for  each  of  the  100  pixels  is 
calculated  and  printed  out.  A header  provides  relevant  housekeeping 
data  including  filter  positions,  th.^  ^ype  of  light  source,  the  actual 
pixel  position  being  illuminated  and  the  number  of  frames  included 
in  the  averaging. 

Summary:  After  all  100  pixels  have  been  illuminated,  a summary  sheet 
is  printed.  This  sheet  contains  the  average  count,  minus  dark  count, 
for  each  pixel  while  it  was  illuminated. 

8.2.2  Honeywell  Spatial  Radiometer 

Pre-flight  tests  have  been  conducted  at  AFGL  and  Honeywell  Radiation 
Center  to  determine  the  response  of  each  detector  as  a function  of  source 
radiance,  angle,  and  wavelength.  The  original  data  orocessing  plan  is 
shown  in  Figures  10  and  11^^^.  As  a result  of  changes  in  instrument 
design  it  has  been  possible  to  avoid  extensive  data  processing  in  the 
spectral  response  (V  vs  X)  calibrations.  The  irradiance  calibrations 
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FULL  FIELD  POINT  SOURCE 


r 


V = Detector  Output 

a = Source  Azimuth  or  Mirror  Scan  Position  (DEG) 

B = Source  Elevation  (DEG) 

T = Source  Temperature  (°K) 

f(V)=  tpd  = Ipd  = Normalized  Transfer  Function  (NTF) 

■-N  S 

WHERE:  E = Irradiance  (W/CM^) 

L = Radiance  (W/CM^  /SR) 

E^j  = E at  10  Volts 

L|^  = L at  10  Volts 

E =E/c^^^  = Normalized  Irradiance 

F(a,S)  = ""  Angular  Response 

= Direction  of  Peak  Angular  Response 

^ = ff  (a,B)  dil  = Effective  Viewing  Solid  Angle 
S(A)  = [j(V;^)/f(V^^)  J /[E'X/E^  j = Spectral  Response 

WHERE:  E^  = Input  Irradiance  at  Wavelength  X 

X^  = Wavelength  of  Peak  Spectral  Response 

y\_  = — dX  = Spectral  Response  Width  (uM) 

L(X„,T) 

WHERE:  L(X,T)  = Blackbody  Spectral  Radiance  (W/CM  /SR/uM) 


Figure  11.  Honeywell  Spatial  Radiometer  Calibrations  Legend 
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(V  vs.  E)  have  also  been  reduced  in  scope  for  the  first  flight.  Hence, 
this  section  will  concentrate  on  the  full-field  and  point  source  (angular 
response)  calibrations.  Much  of  the  software  developed  for  these  would 
be  useful  for  the  other  calibrations. 

8. 2. 2.1  Full-Field  Calibrations  - Program  FLFIELD 

Program  FLFIELD  reads  an  input  data  tape  and  determines  the  full-field 
response  in  volts  for  both  backward  and  forward  sweep  directions  (ramp 
up  and  ramp  down)  for  all  40  detectors.  Initially  it  was  planned  to 
use  a 7-track  digital  tape  prepared  by  AFGL's  decommutation  center 
from  analog  tapes  recorded  during  the  experiments;  however,  problems 
associated  with  saturation  of  the  more  sensitive  detectors  at  higher 
radiances  forced  use  of  the  Honeywell  9-track  digital  t.ipes , also 
recorded  during  the  experiments.  Tables  of  temperatures  and  radiances 
correlated  to  experiment  numbers  recorded  on  these  tapes  were  supplied 
by  experimenters  to  be  input  to  FLFIELD, 

Results  are  printed  and  plotted  for  convenient  review  by  the  experimenters 
They  are  also  written  to  disk  for  permament  file  and  archive  storage. 

A feature  of  the  program  allows  direct  access  to  these  permanent  files  in 
a lieu  of  tapes  if  different  plotting  parameters  (such  as  scale  changes) 
are  requested. 

8. 2. 2. 2 Angular  Calibrations  - FOVPEAK,  FOVPLOT 

These  programs  process  input  data  from  7-track  taoes  produced  by  AFGL 
Decommutation  Center  from  analog  raw  data  tapes.  Saturation  problems 


1 

! 

associated  with  the  full-field  calibrations  did  not  occur  here.  Time 
jumps  in  the  data  are  used  to  separate  different  experiments  each  of 
which  is  characterized  by  a value  of  the  elevation  6 as  supplied  in  a 
table  to  be  input. 

Program  FOVPEAK  scans  the  data  and  produces  an  output  file  containing 
the  (oIq.Bq)  and  value  (in  volts)  of  the  peak  response  for  each  detector. 

This  is  then  used  by  program  FOVPLOT,  together  with  the  full -field 
results,  to  map  out  F(a,3).  Plots  of  F{aQ,3)  and  F(a,3Q)  are  produced. 

For  the  first  flight  these  are  to  be  analyzed  by  the  experimenter 
to  estimate  the  field  of  view  solid  angle  Q.  It  is  expected  that  this 
quantity  will  be  deduced  more  rigorously  in  future  flights. 

8.3  File  Design  for  Flight  Data  Base 

The  final  flight  data  base  design  was  conceived  to  meet  a number  of  objectives: 

1)  All  information  expected  to  be  required  for  extended  study  of  the 
performance  of  any  instrument  will  be  in  one  file.  Thus,  for  each  instrument, 
pertinent  Link  1,  Link  2,  Link  4,  trajectory,  attitude,  and  reduced  data 
will  be  available  simultaneously  so  as  to  facilitiate  further  processing  and 
analyses. 

2)  Files  will  be  separate  per  instrument.  However,  the  slower  data  rate 
for  the  UV  instruments  allows  for  combining  all  4 instruments  on  one  file. 

For  convenience  in  processing,  a file  should  readily  fit  on  one  magnetic 
tape. 


3)  Since  the  significant  resolution  of  the  bulk  of  the  data  is  limited  to 
1 bit  in  14,  the  final  data  base  will  use  packing  of  CDC  60-bit  words  with 
scaling  and  offset  for  efficient  utilization  of  storage. 


i 

i 

4)  Since  many  variables  are  required  much  less  frequently  than  the  basic 

1000  Hz  frame  rate,  these  variables  will  be  stored  commutated.  {; 

f . 

i , 

i ' 

5)  File  design  for  different  instruments  will  be  as  similar  as  possible 
so  that  common  generating  and  user  routines  can  be  used  from  instrument  to 
instrument. 

Tables  1 - 5 on  the  following  pages  show  the  file  format  that  is  proposed 
per  instrument.  Each  word  is  15  bits  long,  and  the  range  dictates  the  scale 
and  offset  that  will  apply.  The  telemetered  data  are  described  in  Reference  4. 

Trajectory  and  attitude  data  are  generated  by  SUA  and  may  be  obtained  from 
separate  files  or  stored  commutated.  Raw  signals  will  be  accessible  for 
checking  purposes,  but  the  final  analyses  will  be  based  on  the  smoothed  back- 
ground signals  (SM  BKGD)  and  the  calibrated  and  background  removed  signals 
(CABR).  The  file  contents  are  subject  to  considerable  change,  but  these 
specifications  have  been  reviewed  with  the  researchers  and  present  the  scope 
of  the  problem. 

Table  6 quantifies  the  size  of  each  data  base.  For  each  IR  instrument  a 
reduction  in  the  output  subframe  rate  is  proposed  either  by  averaging  or 
taking  one  of  two  samples.  Assuming  15  bit  values  or  4 per  CDC  word,  the 
worst  case  of  9.45  x 10^  values  with  the  Honeywell  Spatial  Radiometer  can 
be  stored  on  one  9-track  1600  bpi  tape. 
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NSNP  ENVIRONMENTAL  DATA  BASE 

i/r/ft 

OUTPUT  FILES 

photometers  «N0  SPtCTROW£TER  FILE 


M093  2 IS  A CQHWUTtTOR  <iUll  X It  WHICH  CQHTtMS  t9Y  SFIgM  TJMgji; 
TAUGEf  HIC»G,  2-li;  SENSOR  WPGf  12-20?  SP  CABR,  21-56?  SP  8KG0, 

52-92?  HSP  HKPG.  93-99?  SP  MRPG.  100. 

HORO  3 IS  A CONNUTATOR  HHICi?  CONTAINS!  (OPTIONAL? 

TRAJ  L ATT,  _1 -961  SP  RAM.  97-62?  MSP-1  HKPGj  63-69}  HSP-3  HRPG, 


90-96?  SP 

HKPGf  97»100. 

MONO  SFIO 

SYMBOL  OESCRIPTION 

range 

1 

SFIO  (100  SUBFRAMES/FRAME) 

1/100 

2 1 

TIME  FROM  launch  (SEC) 

0/601 

T_<.  65  SEC 1.^65  . SEC 


I Z TIMER  EVENTS TlNgl^  E»6WTS 


2 

3 

ENGINE  OXIDIZER  PRESS 

0/5 

(VOLTS) 

2 

4 

ENGINE  CHAMBER  PRESS 

0/5 

(VOLTS) 

2 

5 

ENGINE  FUEL  PRESS 

0/5 

(VOLTS) 

2 6 timer  EVENTS 


2 

2 

2 

2 

7 

8 

ENGINE  GN?  PRESS  (VOLTS) 

ENGINE  REGUL  PRESS  iVOLTS) 

0/5 

0/5 

0/5 

0/5 

9 

10 

ENGINE  FUEL  TEMP  (VOLTS) 

ENGINE  OXIDIZER  TEMP  (VOLTS) 

2 

li 

ENGINE  6N2  TEMP  (VOLTS) 

0/5 

2 

12 

DRO  FUNCTION  MON  (VOLTS) 

0/5 

2 

13 

OLD  FUNCTION  MON  (VOLTS) 

0/5 

2 

14 

87A  BATTERY  (VOLTS) 

0/5 

-2_ 

1 5 

B9A  BATTERY  (VOLTS) 

0/5 

2 

16 

K5  LOGIC  RELAY  (VOLTS) 

0/5 

2 

17 

<6  LOGIC  RELAY  (VOLTS) 

0/5 

2 

lo 

K9  LOGIC  RELAY  (VOLTS) 

U /5 

2 

19 

2 

20 

2 _ 

21- 

SP  CABR  (PNOTnNS/CM2/SEC/A) 

56 

2 

57- 

SF  SM  BKGD  <PH0T0NS/CM2/SEC/A) 

92 

2 



HSP  FILTER  * ..  _ . . 

Table  1.  UV  Instruments  Data  Base 


J 


142 


T 


BtSl  AVAILABLE  CO?V 


28 


3 

29- 

TA^lGEJ  ACC£L  tSCSl_  (KM/5EC2> 

-.01/ 

31 

• Cl 

3 

32- 

tJLERIAN  ANGLES  OEF  ECI  TO  SCS  10E6» 

u/360 

3 

34 

35- 

eulepian  Angles  oef  ven  to  scs  (Degi 

0/36C 

37 

3 

3tt- 

TlRiJST  direction  IN  SCS  (D£i,) 

0/36  J 

39 

J 

<*0 

ROLL  angle  3F  TARGET  <OEG) 

u 

3 

41- 

TARGET  NOZZLE  POS  IN  SENSOR  PLANE  (OEG) 

-10/  10 

82 

3 

43 

angle  between  LOS  i HARO  EARTH  TANGENT 

(OEG) 

-90/90 

3 

_>t‘i 

TANGENT  HT  OF  LOS  FROH  NFARgST  EARTH  SORF. 

0/3  JC 

(KH) 

3 

DISTANCE  Along  los  to  pt  nearest  earth 

SURF  . 

0/2500 

3 

46 

target  THRUST-VEL  angle  (0EG> 

0/183 

3 

ht- 

82 

SP  RAW  (CQLINTSI  _ . 

1“^  BITS 

3 

83 

. HSP-1  filter  » 

3 

34 

WSP-1  TENP  A (rfOLTS) 

0/5 

3 

85 

f1SP-l  TENP  B IVOLTSI 

C/5 

? 

86 

M5P-1  TEMP  C (VOLTS) 

0/5 

3 

32 

“ISP-1  H.v. COP'D  (VOLTS) 

C/5 

3 

dd 

MSP-1  LAMP  volt  (VOLTS) 

0/5 

3 

3i . . __ 

MSP-1  CORONA  OET  (VOLTS) 

u/5 

Table  1.  UV  Instruments  Data  Base  (continued) 
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2 

94 

HSP  TEMP  A (VOLTS) 

0/5 

2 

95 

MSP  TEMP  a IVQLIS) 

0/5  ' 

2 

95 

HSP  TEMP  C (VOlTS) 

C/5  j 

2 

9Z 

MSP  M..VAlLm'Il.  (VQLISJ 

0/5  1 

2 

98 

HSP  LAMP  volt  (VOLTS) 

0/5  1 

2 

99 

HSP  CORONA  volt  (VOLTS) 

0/5  , 

2 

130 

SP  TEMP  A (VOLTS) 

0/5  j 

TRA  JECTOT  Y 

1 ATTIIUOE  (WORD  3,  SFIO  1-46)  - OPT  I ONAL 

3 

i 

- RI  ASCER  OF  GREENWICH  lfl££l 

0/360 

3 

2-4 

GEOC  POS  COMPONENTSCSENSOR) 

(KH) 

2500/ 

65  u 0 

3 

5-7 

GEOC  VEL  COHPONENTSC SENSOR) 

(KM/SEC) 

-2/2 

3 

8-lQ 

GEOC  ACCtL.  JlftaPOREMISISENSOR)  (RH/SEC2) 

-.5/. 5 

3 

11 

GEOO  alt  above  SUB-EARTH  PT (SENSOR)  (KH) 

U/30C 

3 

12- 

IU.QC^POS_j:QJlPJlNtNTilT  A RGE  T ) 

(KH) 

2500/ 

14 

6 5 0 0 

3 

15- 

GEOC  VEL  COMPONENTS(TARGET) 

(KM/SEC) 

-2/2 

1 

17 

. la.- 

G/nC  ACCEL  COMPONENTS)  lARGf 

0 (KM/5EC2) 

-.5/. 5 

20 

3 

,21 . 

GRQfl  AO  _AflQlt£_SU3-£ARTH  PT  (TARGET)  (KM)  , 

0/30C 

3 

22 

TARGET  RANGE  (SCS)  ( KN) 

0/15 

3 

23- 

-15/15 

25 

3 

26-  , 

_ TASiiET  VRR  iSCSIs  (Rp/ SEC  J 

-.5/. 5 

BEST  AVAILABLE  COPY 


5 

9aI 

HSP-3  FllltP  • 

i 

91 

MSP-3  TtiF  A (VOLTSI 

t./5 

t 

9C 

rtSP-I  TtHP  a (V0LIS» 

G/6 

i 

9 3 

fiSP-3  TtHF  C IVOLTS) 

n/5 

i 

ISP'S  M.V.COhM  (I»0LTS» 

fi/5 

3 

95 

HSP-3  t»NP  VCtI  (WOLtS) 

1/5 

i 

9b 

•ISP-3  CORONA  VOLI  (VOlTSI 

0/5 

3 

9/ 

SP  TfHP  3 <»(,LIS) 

«/5 

3 

9a 

SP  H.v.rop'o  (Voirs» 

1/5 

3 

93 

SP  LA. IP  volt  (VOLTS) 

-/5 

3 

(* 

S 

100 

SP  C&PONA  OiTTTO-T  (VOLTS) 

HSP  CAjR  (P(CT0NS/C.3?/StC) 

HSP  S I iivol)  ( PHOTONS/CH?/SEC) 

0/5 

7 

6 

3SP  PAW  (COUNTS) 

MSP-1  CA3P  (f  HOTOUS/CM^/SEC) 

MSP-t  sj*  fijc&n  (photons/cn?/sec) 

15  ei TS 

9 

10 

11 

•ISP-1  PAH  (COUNTS) 

NSP-J  CAPf  (PH0T0'(S/CMC/SFC) 

USD-’  SN  aver  (PHOT  JNj/CM?/StC) 

15  BITS 

l^ 

NSP-S  kAH  (C(<UNTS) 

15  BITS 

Table  1.  UV  Instruments  Data  Base  (continued) 
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FULL'FIELO  radiometer  <NR-«)  file 


BEST  AVAILABLE  COPY 


WORD  2 IS  * commutator  (SO  X II  NHXCH  CONTAINS  IBT  SFIOI I TINE, It 
JAR&ET  HXPG.  2-1  It  SENSOR  HKPC.  12-2AI  FFR  MKPC.  2i-2Tt  UNASS.  ZArii, 
NORO  3 IS  A COMMUTATOR  (SO  X 1)  HHICH  CONTAINSI  (OPTIONAL)  TRAJ  A 
ATT.  UNXSR.  hT-MII. 


^MQBtL  SFTD  STMaOL  OESCRIPTIOW 


1 

(50  StBFRAMES/FRANEI 

2 

1 

time  from  launch 

(SEC) 

T < 65  SEC 

T > 65  SEC 

2 

2_ 

2 

TIMER  EVENTS 

TIMER  EVENTS 

2 

4 

( VOL1  SI 

ENGINE  chamber  PRESS 

__2. 

(VOLTS) 

_ ENGINE  FUEL  PRESS  

(VOLTS) 


rank 

i/sg 

0/600 


0/5_ 

0/5 

0/5 


2 

z. 

2 

1. 
2 

2. 


6 TIMER  EVENTS 

Z ENGINE  GN2  PRESS  (VOLTS) 

a ENGINE  REGUL  PRESS  (VOLTS) 


0/5 

0/5 

0/5 

0/5 

0/5 


Z 12 QRO  FUNCTION  HQn  (VOLTS) 

2 13  DLO  FUNCTION  MON  (VOLTS) 

2 15'  ~B10A  3ATTERY  (VOLTS) 

2 16 C2  CAP  MONITOR  (VOLTS) 

2 17 

2 10 KIN  LOGIC  RELAY  (VOLTS) 

2 19 

2 21 

2 21  OT-1  TEMP  (VOLTS) 

2 _-22_ PT-2  TEMP  (VOLTS) 

2 23  OT-C  TEMP  (VOLTS) 

Z Zh PT-H  TEMP  (VOLTS) 

2 25  TB  TEMP  (VOLTS) 

2 Bfl  temp  (VOLTS) 

2 27  FP  TEMP  (VOLTS) 

2 24^ . . 

50 


0/5 

0/5 

0/5 

0/5 

0/5 

0/5 


0/10 
0/10 
0/10 
0/10 
0/10 
0 /lO 
0/10 


TRAJECTORY  t ATTITUK  (H0RD3,  SFIO  l-<i6) -OPTIONAL  ' 

3 1 RT  4SCEN  OF  GREENWICH  (DEG) 

3 2-N  GLOC  POS  COMPONENTS! SENSOR)  (KM) 

3 5-7  GEOC  VEL  COMPONENT  S (SENSO  R)  (KN/SEC) 

i fl-ia GEOC  ACCEL  COMPONENTS ( SENSOR)  (KM/SEC2) 

3 11  GEOD  ALT  ABOVE  SUB-EARTH  PT(SENSORI  (KNI 

3 12-  GtQC  POS  CDIIPONENTS(TARGET)  (AM) 

14 

3 15-  GEOC  ViX  CQMPQNlNTSI TARGE T)  (KM/SEC) 


u /360 
25uD/ 
6500 
-2/2 
-.5/. 5 
0/300 
2500/ 

65  uO 
-2/2 


Table  2. 


Full-Field  Radiometer  (NR-5)  Data  Base 
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best  available  copy 


IT 


J 

3 

3 

3 

14- 

20 

21 

22 

21- 

&£oc  Accet  cohpon:nts<t*rgei)  r>ci</sEC2» 

GEQD  Sy8-£»RTM  PTdARGETI  (<M) 

target  range  (SCSI  (KM) 
target  POS  (SCSI  (KM) 

-.5/. 5 

L / JOO 
u 'IS 
-15/15 

3 

2S 

2&- 

TAI^QET  vet  (SCSI  (KM/ SEC) 

-.5/*5 

3 

?0 

29- 

TARGET  ACCEE  (StS.)  (KM/SEC2J 

-.01/ 

3 

31 

32- 

EJLEKIAN  angles  bee  ECI  to  SCS  (DEG) 

• 01 
0/36C 

3 

34 

35- 

eulerian  angles  oef  ven  to  scs  (Oesi 

u/360 

3 

37 

36- 

THRUST  DIRECTION  IN  SCS  (OEG) 

fl/3bC 

3 

39 

43 

ROLL  ANGLE  OF  TARGET  (OEG) 

b/36 

3 

41- 

TARGET  NOZZLE  POS  IN  SENSOR  PLANE  (OEG) 

-10/  0 

3 

•tZ 

43 

~ "angle  0ETirEEN~LOS  L HARO  EARTH  TANGENT  (OEG) 

-90/90 

3 

<i<i  . 

TANGENT  HT  JQF  _t.OS  PROH  NEAREST  EARTH  SURF. 

u/3ao 

3 

<•5 

(KHI 

DISTANCE  ALONG  LOS  TO  PT  NEAREST  EARTH  SURF. 

0/2500 

3 

4b 

<KH) 

TARGET  THRUST-VtL  ANGLE  (OEGI 

G/iao 

3 

U 

47- 

50 

CARR  (SIGNAL  1 LO)  (M/CM2/SRI 

£-12/ 

S 

_ _ CAOR  ISIGNAL  1 HII  (H/CM2/SR) 

E-4 

£-12/ 

6 

_ CABR  (SI&NA_t  2 MI)  (N/CM2/SR) 

£-4 

£-12/ 

7 

GA8R  (SIGNAL  2 LO)  (H/CM2/SR» 

£ -4 

£-12/ 

ft 

SN  BKGD  (SIGNAL  1 LO)  (W/CM2/SR) 

£-4 

t-lt  / 

9 

SM  BKGD  (SIGKL  1 HI)  (M/CM2/SR) 

E-l  / 

1jL_ 

SM  BFGD  (EIGfAL  2 HI)  (H/CH2/SR) 

E-  2/ 

11 

SM  6KQ0  (SIGNAL  2 UP)  (H/CM2/5R) 

f-b 
£-  2/ 

12 

KAH-OEHOO  L FLTiO  (SIGNAL  1 LO)  (VOLTS) 

J/lQ 

13 

RAh-OEMOD  K FLTRD  (SIGNAL  1 HI)  (VOLTS) 

G/IQ 

Hi 

RAH-DEMOD  i FLTRO  (SIGNAL  2 HI)  OTOLTS) 

0/10 

15 

RAW-OEHOO  i ELTRO  (SIGNAL  2 LO)  (VOLTS) 

0/10 

16 

ram  (SIGNAL  1 LO)  (VOLTS) 

-IC/i 0 

17 

RAM  (SIGNAL  1 HI)  (VOLTS) 

-lt/1 j 

i« 

RAR  (REFERENCE)  (tfOLTS) 

- 1 0 / 1 u 

19 

ram  (SIGNAL  2 HI)  (VOLTS) 

-10/10 

24.. 

RAM  (SIGNAL  2 LO)  (VOLTS) 

-10/1  0 

Table  2.  Full-Field  Radiometer  (NR-5)  Data  Base  (continued) 
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— s 

best  AVWUBlt  COPt 

-VC  SPECTROHCIf >«  SS-6  FILE 


MORO  2 IS  A COHHUIAIOR  (lOfc  X 1)  WHICH  CONTAINS  (BY  SFIOll  TINE,!! 
TARtlll  HILP6.  2-lii  SENSOR  HKPt*  12-2fl;  SCAN  (>  21!  NS-6  HKPC,  22-28! 
LINASSi  29-50:  «OPTIONAL)  tRAJ  ( ATI,  51-96!  UNASS,  9T-1G0. 


HORO 

1 

Z 

SFID 

1 

SYMBOL  OESCRIPIICN 

~ SFIO  (100  subfrahes/franei 

IIN£  from  t.»UNCi(  (S£:> 

T < 65  SEC  I > t5  SEC 

RANGE 

l/iOQ 

u/bao 

2 

2 

TIMER  £!I£NTS  . J INER  EVENTS 

2 

3 

ENGINE  OXIDIZER  PRESS 
( VOi  T S » 

0/5 

2 

w 

ENGINE  CHAMBER  PRESS 
•VOLTS) 

3/5 

2 

5 

ENGINE  FUFL  PRESS 

(VOLTS) 

0/5 

2 

TIMER  EVENTS 

2 

7 

LNGINE  GN2  PF-ESS  (^OLTS) 

0/S 

p 

A 

ENGINE  REGUL  PRESS  (VOLTS) 

c/5 

2 

9 

ENGINE  FUEL  TEMP  (VOLTS) 

3 /S 

2 

13 

_ iNGINt.  OXIOIZtR  T£NP  (VOLTS) 

C/S 

2 

1! 

ENGINE  GN2  TIMP  (VOLTS) 

0/s 

2 

12 

OPi)  FUNCTION  ION  (VOLTS) 

a/s 

2 

13 

Olu  function  MON  (VOLTS) 

s/s 

2 

14 

99ft  RftTTEPy  (VOLTS) 

0 /s 

2 

IS 

810  A 3AT  T£RY  (VQLTS) 

0/s 

2 

C3  CAD  monitor  (VOLTS) 

0/s 

2 

1^ 

•7 

p 

16 

19 

KIO  -OGIC  relay  (VOLTS) 

0/5 

4. 

20 

2 

21 

SCAN  # 

1/12U0 

2 

22 

PT  TCMP  (VOLTS) 

0/10 

2 

23 

OT-C  TLMP  (VOLTS) 

0/10 

2 

24 

OT-W  TEMP  (VOLTS) 

D/lO 

2 

25 

19  TlMP  (VOLTS) 

a/10 

33  Tfhp  (VOLTS) 

0/1  u 

2 

2 

2 

Z7 

2« 

29- 

FP  TEMP  (VOLTS) 

0/10 

53 


TRA JtCTORY 

( ATTITUOt 

. (WO»sO  2,  SFIO  51-96)- 

OPT lONAL 

* 

51 

PT  ASCEN 

OF  GfitENMICH  (0E5) 

o/  360 

2 

S2- 

&£  OC 

POS 

rOHFONc.NTS(SENSOP) 

(<H> 

2^1, 0 / 

54 

oS  0 G 

2 

55- 

GE  )C 

VEL 

COMPONENTS (St  NSCP) 

(K H/StC) 

-2/2 

n 

2 

3 8- 
63 
o 1 

GtOC' 

aCCc.L  COMPONENTS  (SENSOR)  (<H/SEC2) 

-•5/. 5 

2 

GlOJ 

Alt 

ABOVE  SU3-EAPTH  PT  (SENSOR)  (XM) 

a/3Jj 

2 

52- 

GEOC 

POS 

CONPONENTSITARGET) 

(KM) 

25a3/ 

Table  3.  CVF  Spectrometer  NS-6  Data  Base 
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BEST  AVAllAS'i  COPY 


(i£OC  VEl  components «T*IICET|  «ICM/SCC) 
iiEDC  „aqq£i,  _cohpqncnts(T»r&£T»  (KN/sec?> 
ABOVE  sub-earth  PHTARtET)  (KHI 


target  range  (scs)  (khi 
TARfttT  M5  «CSJ  (<NJ 

target  VEE  (SCSI  (<M/SEC» 

target  ACCEL  (SCS)  (KH/SECZ) 

EULERIAN  ANGLES  DEE  £CI  TO  SCS  (DE6I 

EULLRIAN  ANGLES  OEE  VEN  TO  SCS  (OEG» 

THRUST  DIRECT  ION  IN  SCS  (PEG) 

ROLL  ANGLE  OF  TARGET  (OEGT 

TARGET  NOZZLE  POS  IN  SENSOR  PLANE  (OEGt 

anglI  between  LOS  s harY)  tarth  tangent  (OcG) 
tangent  HT  JF  LOS  FROM  NEAREST  EARTH  SURF. 

( <H) 

DISTANCE  ALCng  LOS  TO  PT  NEAREST  EARTH  SURF. 

(KH) 

TARGET  THRUST-VEL  ANGLE  (DEG) 


X OF  SCA^l  _ 

wavelength  (HU) 

CABR  (SIGNAL-LO)  ( WZCM2/SR/Hg) 

CABR  (SIGNAL-HI)  (W/CMZZSR/'NU) 


SN  SKGCL.(-SIGN4L-L0)  (W/CM2/SR/HU) 

SM  BKGO  _(_SI&NAL-HI)  ( W/CHZ/SR/NU) 

RAW  (SIGNAL-LO)  (VOLTS) 

RAW  (SIGNAL-HI)  (VOLTS) 

posnroN  reference  (voltsi 


Table  3.  CVF  Spectrometer  NS-6  Data  Base  (continued) 


Gsao 

-2/Z 


• /J#0 


S/IS 

-15/lS 


C/J63 

-10/10 


-Bu/9o 

0/3)0 


0/100 

Z/6 

E-IZ/ 

E-w 

i-lZ/ 


E-G 

£-12/ 

c-4 

E-12/ 

E-4 

-10/10 
- 10/10 
-10/10 
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C^r  SPfCmOrftUrt  HS-i  riLk 


w mUBLl  copy 


NU^n  ^ lb  A COMHUlAKVt  (lUu  * 1)  MHICH  CONlAINS  (IT  SFin)!  TIHF.i: 
lARGtT  HKP&t  2-li;  StNSU^<  HKPj  , l^-^0:  SCAN  •,  2i;  HS-J  rtKPG,  ZJ-Z#; 
UNASS,  Z9-Si.;  (OPtlONAL)  t>-AJ  1 All,  Sj-96:  UNASS,  97-100. 


UOtiO 

SFIO 

SYM90L  OESCNIPIIGN 

RANGE 

1 

SFIO  (100  SUOF-(AHtS/FRAME» 

• /iJc 

1 

TIME  FROM  launch  (SEC) 

5/60C 

I < bS  SEC  T ► 65  SEC 

a 

2 

TIMEP  EVlMS  timer  FV:lHTS 

t 

3 

engine  OTIOlZtR  PPcSS 

j /5 

(vni IS) 

j 

ENGINE  CHEMndR  PR'^SS 

J/5 

(TOE  ESI 

S 

EIIGlHt  FJ-L  “B'.SS 

u/5 

(VJL  IS) 

z 

6 

TIMER  '•vents 

? 

T 

engine  GNC  press  (VOLTS) 

j/5 

2 

6 

^NGINu  REbUL  PRESS  (VOLIS) 

0 /5 

? 

9 

ENGINE  FUtL  temp  (VOLTS) 

G/5 

2 

ill 

engine  OTIDIZER  temp  (VOLTS) 

0/5 

2 

1 1 

ENGINE  &N3  Tl  MO  (VOLTS) 

,a/0 

? 

12 

ORO  FUNCTION  MON  (VOLTS) 

0 /5 

7 

1 J 

CLP  FUNCTION  MON  (VOLTS) 

0/5 

2 

1 4 

99A  aATTEPy  ( VOLTS) 

0/5 

2 

15 

dlOA  MATTERY  (VOLIS) 

0/5 

2 

1 o 

r,4  CAP  MONITOR  (VOLTS) 

J /5 

2 

17 

<10  logic  relay  (VOLTS) 

0/5 

2 

IH 

<14  logic  relay  (VOLTS) 

0/5 

2 

19 

2 

20 

2 

21 

scan* 

1/I2o0 

2 

22 

pre-temf  (vnL[S) 

0 /ic 

2 

25 

J-T  TEMP  (VOLTS) 

3/10 

7 

24 

CF-C  temp  (VOLTS) 

3/lJ 

2 

25 

CF-N  TEMD  (VOLTS) 

3 /l  u 

2 

26 

BAF  TEMP  (VOLTS) 

0/10 

2 

27 

OOT  TSmP  (VOLTS) 

0/1^ 

2 

25 

2 

29- 

S3 


T<AjFGI1'<T  ( ATT1TU9-:  (HO'.J  ?,  SFIU  SI -9b  ) -uPT  1 JNAL 


2 

51 

<l  AS 

CFN 

t.'f  (135  f> 

/ i t,'  « 

5 

3 ’ - 
5 u 

Gc  00 

POb 

roMPnNLNTS(ifN53Pi 

(Kl) 

25o  1 / 
bS 

7 

^9- 

. 7 

GcUC 

rO  -lE  O'i^MTS  ( N:>nC'j 

(<  •VSi.C) 

-2/2 

J ! 

'>i- 

5C 

ol 

G-  OC 

acc* 

:l  *11  s ( 

p)  K<'^/i>'ZZ\ 

- . 5 / » 

2 

m 

ALT 

fl90V_  5U9-':a=TH  FI 

(b* 

./■'wl 

? 

•j2- 

0 0 

pn:> 

CP'IK  9N_NT  S ( I fl  1 ) 

(<  1) 

J5j3/ 

Table  4.  CVF  Spectrometer  HS-3  Data  Base 
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BtSI  AVKIE  CC?lf 


r 


i 

64 

6S- 

IjlOC  V£W  COMPONENTS!  TAR&ET)  !KM/S£C> 

bSfiO 

-2/2 

1 

2 

67 

6«- 

70 

71 

Otoe  *C£tk  COMPOM£NTS<T*RGET)  <KH/SEC2) 

GEOO  ALT  ABOVE  SUB-EARTH  PT  (TARGET)  (KNI 

-.5/.5 

ai/3i}0 

2 

2 

2 

72 

73- 

75 

7b- 

TARGET  range  (SCS)  (KM) 

TARGET  POS  (SCS)  (KM) 

TA^fcT  VEL  (SgS)  (KM/SEC) 

0/19 

-19/15 

-. 9/. 5 

78 

79- 

TARGET  ACCcL  (SCS)  (KM/SEC2> 

-.Ml/ 

2 

81 

82- 

EUUERIAN  angles  DEF  iCl  TO  SCS  (OEG) 

. 01 
0/160 

2 

84 

85- 

EULERIAN  angles  DEF  VEN  TO  SCS  (OEG) 

m/360 

2_ 

87 

- e&-  

thrust  direction  in  SCS  (OEG) 

0/36C 

2^ 

69 

90 

ROLL,  angle  (}F  target  (OEG) 

0/360 

2 

91- 

92 

TARGET  NOZZLE  POS  IN  SENSOR  PLANE  (DEC) 

-10/10 

2 

2 

93 

94 

ANGLE  between  LOS  E HARD  EARTH  TANGENT  (OEG) 
TANGENT  HT  OF  LOS  FROM  NEAREST  EARTH  SURF. 

-90/90 

0/  JOO 

2 

2 

99 

96 

(KH) 

distance  along  LOS  TO  PT  NEAREST  EARTH  SURF. 

(KHI 

TARGET  THRUST-VEL  ANGLE  <0£GI 

0/2500 

0/160 

2 

9^- 

3 

100 

V.  OF  scan 

Q/lOO 

4 

5 

WAVELENGTH  (HU) 

CABR  (SIGNAL-LC)  ( W/CMZ/SR/HU) 

2/6 

£-12/ 

6 

CA8R  (SIGNAL-HI)  ( WZCMZ/SR/MU) 

t- 

£-12/ 

2 

Si  3KGD  (SIGNAL-LO)  < W/CH2/SR/HJ) 

E-4 

£-12/ 

8 

SM  BtCGO  (SIGNAL-HI)  (W/CM2/SR/HU) 

t-4 

E-12/ 

._2_ 

RAW  (SIGNAL-LO)  (VOLTS) 

E-h 

-10/10 

10 

11 

RAW  (SIGNAL-HI)  (VOLTS) 

POSITION  reference  (VOLTS) 

- 10/10 
■ lO/lO 

Table  4.  CVF  Spectrometer  HS-3  Data  Base  (continued) 
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BEST  AVAIl/iBlE  COPY 

HONtYUFlL  SP»!I«l  8»0I0I*TER  FILE 


M)R3  Z IS  « COnHUTETOR  (ZS  X 11  MHICH  CONTtINS  IBT  SF 1 0>  I TIHEt  II 
TARGE.!  HRPG.  Z-HI  SEHSQR  HRPG.  \Z-Zil  SHEEP  ti  Zll  UNASS i ZZ-Z5. 

MORO  3 IS  A CONHUIATOR  WHICH  CONTAINSI  lOPTIONAL)  TRAJ,  1-ZS. 

HQRP  1.  IS  A CQNWUTATQR  WHICH  CONTAINS!  (QPTlONALi  TRAJ  A ATT,  l-Zj} 

UNASS,  ZZ-Z5. 

HORO  SFIO  SYHBOL  OtSCRIPIION  RANGE 

1 SFID  (Z5  SUBFRANES/FRAfC)  1/ZS 

Z 1 TIME  FROH  LAUNCH  ISEC) il/biO 

I < feS  SEC  T > 65  SEC 


z 

2 

HH£R  EVENTS  TIHER  EVENTS 

z 

3 

ENGINE  OXIDIZER  PRESS 

0/5 

(VOLTS) 

2 

4 

ENGINE  CHAMBER  PRESS 

0/5 

(VOLTS) 

2 

5 

ENGINE  FUEL  PRESS 

0/5 

(VOLTS) 

2. 

h 

timer  events 

Z 

7 

cNGINt.  GNZ  PRESS  (VOLTS) 

0/5 

z 

8 

ENGINE  RiGUL  PRESS  (VOLTS) 

0/5 

z 

9 

ENGINE  FUEL  TENP  (VOLTS) 

0/5 

•> 

IJ 

ENClNt  oxidizer  TEHP  (VOLTS) 

0/5 

11 

ENGINE  GNZ  TEHP  (VOLTS) 

0/5 

2 

12 

ORD  FUNCTION  HON  (VOlTSI 

u/5 

2 

13 

auo  function  «on  (yoLTSJ 

D/5 

? 

14 

811A  3ATTERY  (VOLfSJ 

0/5 

z 

1^ 

Cl  cap  hohIToR  (vouts) 

0/5 

2 

Xt> 

KIO  LOGIC  RELAY  IVOLTS) 

0/5 

l_ 

ir 

Kll  LOGIC  RFLAv  (VOLTS) 

0/5 

2 

t 9 

2 

2 

2J 

2 

21 

SVE.2P  • 

I/bOO 

2 

22- 

25 

I RAJ  L AT  I 

())ORO  i,  SEID  1-ZE  and  WORO  4,  SFJD  1-ZI) -OPT  lONAL 

3 

1 

RT  4SCEN  OF  GREENWICH  (DEG) 

0 /3bu 

i 

Z-R 

GEOC  PUS  COMPONtNTSISENSOR)  (KM) 

25-0/ 

6500  . 

5-T 

GEOC  VEL  COHPONENIS(SENSOR)  (XH/SEC) 

-2/2 

3 

d-13 

&ECC  ACCEL  components (SeNSOR)  (<H/SEC2) 

-.5/,5 

11 

GiOO  alt  above  SIB-EARTH  PKSENSOO)  (KH) 

0/33C 

3 

12- 

GEOC  POS  COMPONENTS(TAPGeT)  (KM) 

2503/ 

14 

6500 

3 

IS- 

GEOC  VtL  COHfONLNTSiTARGLT)  (KH/StC) 

-2/2 

.J.?  _ 

3 

1 3- 

GEOC  ACCEL  COHPONENIS(TARGET)  (XP/SECZ)  ~ 

-.5/. 5 

20 

i 

21 

GEOO  alt  ABOVE  SUB-LARTH  PT  (TARGET)  (KN) 

j/ J JO 

3 

22 

IARjET  range  (SCS)  (RD 

0 /15 

Table  5.  Honeywell  Spatial  Radiometer  Data  Base 
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BtSl  AVAIIABIETOPY 


0 

2 ?- 

lAR&tl  FOb  (SCS)  «KM» 

- 19/  1 9 

4 

1-3 

TAKOtl  WCL  <iC5»  (<1/SEC) 

-.9/.S 

4 

4-b 

TAkUET  ACCtL  (SCS)  (KM/S£CZ» 

-.01/ 

• 01 

4 

7-S 

FlILKRIAN  ANM'e'  IJ£F  ECI  TO  SCS  (OtO) 

U/  360 

4 

10- 

EJLt'^IAN  AJGLts  O'F  VE  N TO  SCa  (DFj» 

0 /3bl 

12 

4 

13- 

THBUSI  DIRECT  ION  IN  SCS  ( OEGt 

w / 3 b 0 

4 

1 4 

IS 

ROLE  ancle  Of  target  tOEG) 

0/  3b0 

4 

16- 

target  NOZZLE  POS  IN  SENSCP  FLANt  (0£G» 

-IJ/iO 

1 7 

4 

19 

ANGLE  OFTWEEfi  LOS  L HARO  EARTH  TANGENT  (TEG) 

-9U/90 

4 

19 

tangent  HT  Of  LOS  FROM  NEAREST  EARTH  SURF. 

u/30C 

(KH  ) 

4 

2J 

DISIANCt  ALONG  LOS  TO  PT  NEAREST  EArTH  SURF. 

u /29u  u 

(<H) 

4 

21 

TARGET  THKIIST-WEL  ANGLE  (OEGI 

u/lOw 

2S 

S 

1 -6 

HK  1 -HSrt 

S 

9-16 

TMPl-TMPS 

17- 

6 HI-<ROP  POS 

C^^^<  FPUH  £Af  H OF  40  DtT^CTOPS 
47-  SM  %JiO  F"UM  =.flCH  OF  4J  niTFCTO-?S 

9*0 

i/-  KAH  F^OM  fAO  OF  4C  OcTFClO^S 

IJf* 


Table  5.  Honeywell  Spatial  Radiometer  Data  Base  (continued) 
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' Section  9.  General  Programs 


Author:  K.  H.  Bhavnani 


9.0  Packing-Unpacking  Programs 

Initiator:  E.  Robinson 

Project  No:  2123  Problem  No:  4937 

Two  general  subroutine  systems  were  developed  for  efficient  packing/unpacking 
and  input/output  of  fixed  length  records.  They  have  been  employed  in  various 
instances  where  large  quantities  of  data  in  repeating  format  must  be  stored 
and  processed.  Subroutine  PACKER  handles  records  consisting  of  an  independent 
variable  (such  as  time)  requiring  a full  CDC  60-bit  word,  and  a fixed  number 
of  NB-bit  words  where  NB  is  a sub-multiple  of  60. 

Subroutine  PACKUP  handles  records  consisting  of  a fixed  number  of  variable  bit 
length  words,  the  only  requirement  on  the  bit  length  being  that  packed  data 
must  fit  60-bit  boundaries.  All  unpacked  data  are  taken  to  be  contained  in 
the  rightmost  bits  of  CDC  words;  therefore  it  is  incumbent  on  the  user  to 
scale  his  data  depending  upon  range  and  resolution  limits  so  as  to  always  re- 
present them  as  non-negative  integers. 

User's  guide  listings  for  the  available  routines  are  provided  on  the  next  two 
pages. 

9.1  Functional  Description 

Packing  and  blocking  procedures  are  used  in  order  to  optimize  for  storage  and 
speed.  BUFFER  IN-OUT  is  used  since  this  is  faster  than  Fortran  binary  read- 
write,  and  to  allow  simpler  information  interchange  with  other  computer  in- 
stallations if  required.  510-word  blocking  is  used  since  this  is  an  accepted 
block  length  for  800  bpi  tapes,  and  to  remain  within  the  512-word  binary  block- 
ing of  SCOPE  3*4  even  if  Fortran  binary  read-writes  were  desired  for  some 
reason . 
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BEST'AVAILABLETOPY 


{ 

( 

I 

( 

( 


c_ 

SUMOUTINE  PACKER  PACKL 1 9.  lOwL C6IC0M.  MAY, 

19ZZ 

C 

c 

SURROIITIMF  PACKER  PROCESSES  FIXED  LENGTH  RECORDS  CONSISTING  OF 

c 

TIME  AND  THE  FIXED  RIGHTMOST  *F-B*  BITS  OF  'LO*  DATA  NOROS.  THESE 
PACKED  INTO  'MH*  COC-HORD  RECORDS.  AND  INPUT-OUTPUT  TO  SEQUENTIAL 

ARE 

ETLFC. 

c 

C 

THE  FIRST  NORD^(EG.  TIME)  IS  STORED  UNCHANGED.  AND  THE  RIGHTMOST 

N8  "I^S 

c 

£_ 

OF  THE  LO  DATA  WORDS  ARE  PACKED  FROM  LEFT  TO  RIGHT  INTO  COC  60-eiT  WORDS. 
WTTM  THF  UNUSED  PART  OF  THE  LAST  CDC  WORD  OF  EACH  RECORD  ZERO-FILLED. 

C 

c 

EACH  RECORD  THUS  CONSISTS  Of  H4  « 1 ^ LO^NB/bO  <ROUNOEO-UP)  COC 
AN  TMTEGRAL  NUMBER  -NR*  OF  THESi  RECORDS  ARE  BLOCKED  INTO  510  COC 

c 

THE  FIRST  HORO  OF  EACH  BLOCK  IS  A CONTROL  WORD  WHICH  DEFINES  NB, 
THF  UNUSED  WORDS  AT  THE  END  OF  EACH  SLOCK  ARE  ALSO  ZERO-FILLED. 

NW  X NO. 

c 

r 

actual  I-O  is  ACCOMPLISHED  BY  BJFFER  IN-OUT  STATEMENTS  USING  000 

'^APIT  Y. 

ALL  ROUTINE  CALLS  ARF  OF  THf  FORM 
C CALL  AYIPEITAPE.  TIMC,  EARRAE.  lN.  LARRAYt  , RETURNS  IE.  F ) 

C WHERE 

C_  ITAPE  IS  THE  LOSICAl  FILE  NUHR^R  EFROH  X TO  10  ALLOWED) 

C lARRAY  ANO  LARRAY  MUST  HAVE  OIMi-><3Ion  lO 
IN  HAY  TAKE  ON  ANY  VALUE  FROM  1 TO  LD 


C E IS  THE  ENO-OF-FILE  RETURN. 


c IS  THE  PARITY  ERROR  RETURN 


THE  list  of  available  ROUTINFS  ROLLOHSI 


call  MUTPtITAPE,  NO,  LD) 

ESTATLISHFS  CONTROL  INFCRHATIO>’  rpi;  wRITInS  OR  READINr.. 

N8  MUST  BE  A PROPER  S'lO  - MUL  T I P. OF  60 

“call  OLTP  (I ’’APE  I 

FQR  READING  PREVIOUSLY  PREP  A RE  D J^^E  S , CONTROL  INFORMATION  MAY  BE 
INITIALIZE!.  FROM  THf  )LOCK  THAT  15  NEXT  IN  POSITION. 


C CALL  PKTPCITAPE,  TI-E,  IARRAY,  LN) 

C TIME  AND  THE  first  LN  WORDS  fr.'m  IARRAY  ARE  PACKED  OUT. 

C 

C CALL  UPTPIITAPE,  TIMF,  IARRFY,  LN,  01,  RETURNStE. FI 

C TIME  AND  THE  FIRST  LN  WORDS  ARr  UNPACKED  INTO  IARRAY.  “ 

C 

C CALL  WRTPdTAPE,  TIME,  IAPO.'y,  lN,  LARRAY) 

C LARRAY  SHOULD  CONTAIN  LN  INTFGi'RS,  EACH  HAVING  A VALUE  FROM  1 TO  LO. 
C THE  CORRESPONQING  IARRAY  LOCATtqNS  ONLY  REPLACE  EXISTING"  R'aCKEO  OA  T J 


CALL  ROTPdTAPE,  TIME,  lARRry,  lN,  LARRAY),  RETURN3(E,F) 

TIME  AND  THE  SELECTED  LN  HOROS  ONLY  ARE  UNPACKED  INTO  I A I P A Y ( 1 -LP ) , 


CALL  OUTPdTAPE.  JIMF^  lARR^rJ  _ _ 

TIME  ANO  NN  PRE-PACKFO  WORDS  ARE  COPIED  OUT  FROM  lAPRAY. 

call  INTP(ITAPE“i  TIME,“iAPPtY,“  0,  0),  RE  TURNS'lE  , F ) 

TIME  AND  ALL  NW  PACKFP  HOROS  A^:  COPIED  AS  IS  INTO  IARRAY. 


call  iSTPflTAPE) 

TACKSPACE  OF  JNL  WECOR’'  ONLY  Wh-n  u.iA)!NG  A FILE. 

CALL  FNTPd  TAPE)"  ‘ 

THF  last  OLOCK  is  WPITTFN  TD  TfPf  FOLLOWED  3Y  AN  LNI.-)F-FILL. 


_C»i^U  SE.lPIlTaPE) _ 

INDICES  Ar'e  resit  by  Onf  to  ALLOW  a"nCW  f"ilF  TC  ACCF  JS  ThFSC  ROuIINi' 
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BEST  AVAILABLE  COPY 


■•c_ 

• C 
■£_ 
- c 

c. 

■ c 

■ C- 

■ c 

• c 

' c_ 

' c 

■ C-_ 

• c 

c 

c 


c 

c 

SL. 
' c 

■ c 

■ c 

c 

' i_ 

■ c 
c_ 
c 
c 

■ c 

■ JL. 
c 

5.. 

• c 


I.  ID-LOGICON. 


MAY,  IHTT 


SU—auTIMg  PACKJP  rWOCesSES  FtYED  LtWGTM  RECORDS  CONSISTING  OF 

THE  iARIABlE  «I6HTH0ST  'MB'  BITS  OF  'LO*  OAT»  NOROS.  THESE  A»C 
PACKEO  INTO  'fu!*  COC-MORD  RECQgQS,  AND  IMPUT-OUTPUT  TO  SEQUENTIAL  FILE' 


TNF  LD  tlATA  WOROS  ARE  PACKED  FPQM  iEflig  RIGHJ  INTO  COC  f.0-BlT  NORn«;. 
MITM  ALL  PACREO  DATA  EYCFPT  THF  LAST  TERNIMAIING  ON  6a-‘lIT  BOUNOARIE^, 
Awfl  TNF  LINUSF3  PART  OF  THE  LAST  CPC  WORD  OF  EACH  RECORJ  7ERO-FILLED. 

AN  INTEGRAL  NUMBER  'NR'  OF  THE'tE  RECORDS  ARE  BLOCKED  INTO  510  COC  MO»OS. 
TNF  FIRST  MQRO  OF  EACH  BLOCK  IS  A CONTROL  NORO  WHICH  OEFIMES  NM  I NP, 

THE  UNUSEO  NOROS  AT  THE  ENO  OF  EACH  BLOCK  ARE  ALSO  ZERO-FILLED. 
tCTUAt  I-O  iS  ACCONPLISHEO  8Y  BUFFER  I?L*5yT  SJATENENTS  USING  OOP  PARITY. 

ALL  ROUTINE  CALLS  ARF  OF  THE  FQBH 

CALL  KYFLdTAPE,  lARRAY,  LN»  RETURNSIE.Ft 

LLHERL—  — — 

ITAPE  IS  the  logical  FILE  NUMBER  ( FF ON  1 TO  10  ALLONEO) 

lARRAY  AND  LAkRAY  MUST  HAVE  OIMENIION  LD  _ _ _ 

LN  MAY  TAKE  ON  ANY  VALUE  FROM  1 TO  LD 

E IS  THE  END-OF-FILE  RETURN,  ^ IS  THE  PARITY  ERROR  RETURN 

THE  LIST  Of  AVAILABLE  POUTINES  ■^OLLOHSI  

CALL  NULLUTAPEj  lARRAY,  LT.'  _ _ . 

ESTABLISHES  CONTROL  INFORMAT  IO*i  FQR  WRITING  OR' REAniNG. 

lARRAY  PROVIDES  THE  VALUE  NB  F QR  EACH  DATA  WORD. 

C^LL  OLFLIITAPE,  lARRAY,  LOt  „ _ _ _ _ _ 

FOR  READING  PREVIOUSLY  RREPARE"’  RILES,  USE  OF  OLFL  RATHER  THAN  NUFL 
VERIFIES  THE  PACKING  INFORMATI-'N  AGAINST  NW  FROM  THE  JJEYT  BLOCK. 

CALL  PKFLdTAPE.  lARRAY,  LSI 

THE  FIRST  LN  WORDS  FROM  lARRAY  ARE  TACKED  OUT, 

CALL  UPFltlTAPET  IAR'rAY,  TniT  RETURN'^rTFl  ^ 

THF  FIRST  LN  WORDS  ARF  UNPACKED  INTO  lARRAY. 

CALL  OUFLdTAPE,  lARRAYl 
C 
C 

c 
C 
c 
c_ 
c 
c 
c 
c 


NM  pre-packed  WORDS  ARF  COPIFO  OUT  FROM  lARRAY. 

CALL  INFL  (ITAPE^  TaRRAY,-  1)1.  RITuRNSIE,  F I 
AWL  NM  PACKED  WORDS  ARE  COPIFO  AS_iS  INTO  lARPAY, 

CALL  BSFLdTAOEl 

BACKSPACE  CF  ONE  PECORO  ONLY  when  REAOInG  1 FILE. 

CALL  enTlTitape) 

THE  LAST  BLOCK  IS  WRITTEN  T"  TApE  rCLLOWEO  jY  AN  FN'i-OE- 
CALL  REFLdTAPFl 


INDICES  ARE  RESET  BY  ONE  T.)  ALLOW  A NEW  FIlF  TO  ACCESS  THESE 


The  first  word  in  each  block  is  the  control  word.  For  PACKER  this  contains 
I zeroes  in  the  first  (leftmost)  12  bits,  block  length  (510)  in  the  second  12 

► bits,  bit  length  NB  in  the  third  12  bits,  number  of  CDC  words  per  packed  re- 

I 

) cord  in  the  fourth  12  bits,  and  the  number  of  packed  records  in  the  block  in 

the  rightmost  12  bits.  This  information  is  sufficient  to  initialize  the 
routines  for  unpacking  a previously  prepared  tape. 

For  PACKUP  the  control  word  information  differs  only  in  that  the  total  bit 
length  for  all  the  variables  is  stored  in  the  third  12-bit  field.  Initial- 
ization before  packing  or  unpacking  requires  that  the  number  of  data  words 
and  the  bit  length  per  word  be  provided.  A feature  of  this  system  is  that  a 
tape  prepared  by  PACKER  can  be  read  by  PACKUP,  the  latter  comprising  more 
general  but  somewhat  (30%)  slower  programs. 

The  routines  allow  various  tape  input-output  operations  after  initialization. 
Pack  (write)  or  unpack  (read)  may  be  performed  in  3 ways: 

1)  an  array  of  the  first  N data  words 

2)  individually  selected  words  from  an  array 

3)  Bulk  transfer  of  the  full  packed  record 
An  ENTP  or  ENFL  operation  is  necessary  at  the  end  of  writing  to  taoe,  to  en- 
sure that  the  last  block  is  output. 

Simultaneous  use  of  either  of  these  routines  is  restricted  to  two  tapes.  RETP 
or  REEL  returns  control  area  storage  for  a new  file  after  a file  has  completed 
using  the  routines.  Logical  tape  numbers  are  restricted  to  1-10.  Finally,  in 
PACKUP  the  maximum  number  of  data  words  is  restricted  to  50. 

9.2  Logical  Procedures  and  Performance 

The  main  performance  consideration  is  that  compared  to  binary  unpacked  READ/ 
WRITE  to  tape,  savings  of  the  order  of  75%  can  be  obtained  in  file  storage  re- 
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qui rements  for  packable  data.  As  shown  below,  data  transfer  times  using  PACKER 
or  PACKIJP  are  also  superior  to  binary  unpacked  Fortran  read-writes. 


Each  of  the  subroutine  calls  requires  only  a few  lines  of  compact  code.  The 
program  is  optimized  in  FORTRAN  IV  by  the  exclusive  use  of  inline  functions 
such  as  AND,  SHIFT,  etc.  Indexing  is  avoided  where  this  is  advantageous.  A 
variety  of  test  runs  were  conducted  leading  to  the  present  implementation  of 
the  programs.  Comparative  execution  times  are  discussed  below,  based  on 
10,000  operations  of  the  type  indicated.  Packing  or  unpacking  times  for  the 
same  type  commands  are  comparable,  and  the  average  values  are  given. 


BUFFER  IN  (OUT)  510  words  3.5  secs 

Binary  READ(WRITE)  510  words  40.0  secs 

All  the  following  handle  time  and  48  data  words. 

Binary  READ  (WRITE)  - No  Packing  8.0  secs 

All  the  following  handle  packed  data  (12  bits). 

PKFL  or  UPFL  7.0  secs 

PKTP  or  UPTP  5.0  secs 

PKTP  (UPTP)  w.  Binary  READ/WRITE  5.5  secs 

instead  of  BUFFER  IN/OUT 


PKTP  (UPTP)  w/o  blocking  11.5  secs 

PKTP  (UPTP)  using  LBYTX,  etc.  19.0  secs 

WRTP  or  RDTP  9.0  secs 

WRTP  (RDTP)  replace  4 of  the  48  data  words  1.7  secs 

OUTP  or  INTP  (direct  transfer  of  the  1.0  secs 

packed  words) 
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9.3  Job  Setup  and  Program  Restrictions 


These  routines  require  about  3K  octal  core  memory 

When  in  frequent  use,  EDITLIB  has  been  used  to  store  the  binaries  on  per- 
manent file. 

Program  halts  due  to  violation  of  restrictions  occur  under  the  following 
conditions; 

STOP  PACKER  ERROR  occurs  if  logical  Tape  number  <1  or  > 10 
or  More  than  2 logical  files  initialized 
or  NB  (bits/word)  is  not  a sub-multiple  of  60. 

STOP  PACKUP  ERROR  occurs  if 

Logical  Tape  number  <1  or  >10 

or  More  than  2 logical  files  initialized 

or  Packed  words  do  not  fit  60-bit  boundaries 

or  Number  of  data  words  >50 

or  Number  CDC  words  calculated  in  OLFL  > control  word  info. 

STOP  BACKSPACE  ERROR  occurs  if  a BUFFER-IN  after  backspace  results  in 
an  end-of-file  or  a parity  error. 

The  BUFFER-  IN/OUT  operations  allow  for  end-of-file  and  parity  error  returns. 
Following  parity  error  on  a read  operation,  the  subroutines  are  set  to  bring 
in  and  process  the  next  block  on  an  ensuing  unpack  call.  However,  a BSTP  or 
BSFL  backspace  call  allows  the  first  record  only  of  the  block  with  the  parity 
error  to  be  unpacked  if  desired. 
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Section  10.  Astronomical  Ephemerides 
Authors:  K.  H.  Bhavnani 
R.  R.  Hayes 


10.0  New  JPL  Planetary  Ephemeris  (DE-96) 


Initiator:  L.  Telford 

Project  No:  4643  Problem  No:  4344 

The  JPL  ephemeris  system  became  unworkable  when  all  copies  of  the  old  JPL  DE19 
tapes  proved  unreadable  for  1977  and  beyond.  JPL  has  gone  through  various 
official  releases  of  Ephemeris  Tape  systems.  The  latest  version  was 
accordingly  obtained  from  JPL^’^.  This  BCD  tape  DE-96  was  converted  to 
binary  using  Fortran  code  provided  on  the  same  tape.  The  new  version  obsoletes 
the  old  JPL  binaries,  is  based  on  updated  planetary  observational  data,  and 
uses  improved  interpolation  techniques  to  compact  data  through  the  year  2000 
on  about  half  a reel  as  compared  to  three  reels  for  the  previous  version. 

Routines  to  process  the  new  tape  were  incorporated  into  System  II  of  the 
JPL  Ephemeris  Processing  System.  This  produces  the  annual  tape  for  Prospect 
Hill  Radio  Observatory.  Results  duplicated  the  earlier  version  for  1976 
to  within  0.5  seconds  of  arc.  A 1977  PHRO  tape  was  then  created.  This 
revision  can  also  be  incorporated  in  System  I - the  general  Ephemeris 
Processing  System,  if  planetary  observations  are  required. 

10.1  Implementation  at  AFGL 

A 7- track  even  parity  magnetic  tape  which  constitutes  the  export  or  users' 
package  for  DE-96  was  requested  from  JPL.  This  encoded  form  is  originally 
prepared  on  a Univac  1108,  but  is  as  machine-independent  as  possible  with 
instructions'^  as  to  how  to  adapt  and  use  the  tape.  There  are  two  files 
on  the  tape: 
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File  1 contains  in  card  image  form  all  the  Fortran  software  needed  for 
conversion  of  the  encoded  ephemeris  file  to  binary,  and  for  use  of  the 
binary  tape.  Twelve  programs  and  subroutines  are  included. 

File  2 contains  the  encoded  ephemeris  in  550  records  of  24000  characters 
each. 

1)  File  1 was  punched  out  to  obtain  a card  deck  for  further  work. 

2)  File  2 was  Buffered-In  in  even  parity  2400  word  blocks  and  written  out 
with  SCOPE  3.4  binary  writes  to  Tape  JPLE96  for  intermediate  storage. 

3)  The  deck  of  programs  was  debugged  for  Fortran  errors  on  the  CDC-6600. 

4)  Program  BCDEPH  and  subroutine  RCI  were  run  using  JPLE96  to  create  the 
compact  binary  tape  JPL  96B  which  is  the  exclusive  data  base  hereafter.  The 
range  of  the  data  is  from  December  16,  1944  to  January  25,  2000. JPL 
achieves  compaction  by  providing  Chebyshev  coefficients  for  polynomials 
which  can  be  interpolated  for  the  ephemeris  of  a particular  body. 

5)  Program  XSHORT  and  subroutines  LH,  SECCAL  and  VIGSEC  may  be  used  to 
create  a short  copy  of  JPL96B  over  a selected  range  of  epochs.  This  feature 
was  tested  but  not  used  in  the  further  development. 

6)  The  main  user  subroutine  is  READE  which  reads  and  interpolates  the  binary 
ephemeris  file.  Double  precision  is  provided  and  has  not  been  eliminated 
though  it  is  not  required  on  the  CDC-6600.  READE  in  turn  requires  sub- 
routines XSTATE,  XINTRP  and  LH.  The  description  of  subroutine  READE  is 
reproduced  from  Reference  2 on  the  following  two  pages. 
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RKADE 

Purpose : Read  and  interpolate  the  binary  ephemerls  file.  Many  users  have  a 

subroutine  by  the  same  name  that  reads  earlier  JPL  ephemerls  files.  This 
routine  can  be  used  in  its  place  without  change.  The  binary  file  must  be 
I on  unit  12. 

I 

■ use : The  user  must  have  the  four  statements 

I 

i DOUBLE  PRECISION  AU,  RE,  TPD,  EMRAT,  TABOUT,  NUT,  JED,  TSEC 

I COMMON/CETBLl/AU,  RE,  TPD,  EMRAT 

C0MM0N/CETBL2/ICW,  ICENT,  IREQ(13) 
i COMMON/CETBL4/TABOUT(6,  12),NUT(4) 

■ in  his  program. 

I 

Then  the  call 

' CALL  READE  (JED,  TSEC,  TERR) 

provides  the  interpolation,  according  to  the  following  rules: 

AU  is  the  scale  factor  that  multiplies  the  output  quantities.  If  AU  = 0, 
all  di.stance  units  are  in  kilometers.  If  AU  = any  other  number,  the  distance 
unit  is  (astronomical  unit)*AU. 

IE  is  not  used. 

TPD  is  the  time  unit.  If  TPD  = 0,  time  units  are  seconds.  If  TPD  = any 
other  number,  time  units  are  days*TPD. 

EMRAT  Is  the  earth-moon  mass  ratio.  If  EMRAT  = 0,  the  value  81.3007  is  assumed. 

ICW  is  a flag  word.  ICW  2 forces  a new  raw  data  read  on  the  next  call  to 
M'iADE.  The  subroutine  sets  ICW  to  2. 

[CENT  1.:  the  center  number  of  the  coordinate  system.  All  vectors  are 
|■'.■f ertinced  to  earth  mean  equator  and  equinox  of  1950.0.  The  bodies  correspond- 
ing ro  center  numbers  are: 

■Number  Body 

1 Mercury 

2 Venus 

3 Earth 

4 Mars 

5 Jupiter 

6 Saturn 

7 Uranus 

8 Neptune 

9 Pluto 

10  Sun 

1 1 Moon 


I 


\ 


I 

I 

I 


I 


i I 


IREQ  Is  what  kind  of  Interpolation  is  requested  for  each  body.  The  meaning 
is : 


IREQ(I)  = 0 , no  interpolation 

= 1,  position  only 

= 2,  position  and  velocity 

For  1=1  through  11,  the  body  designation  is  the  same  as  for  ICENT.  IREQ(12) 
refers  to  the  Earth-Moon  barycenter,  and  IREQ(13)  refers  to  nutations. 

The  calling  sequence  arguments  JED  and  TSEC  specify  the  ephemeris  epoch  at 
which  interpolation  is  wanted.  The  ephemeris  date  is  taken  to  be  JED  + 
TSEC/86400.  Any  combination  of  JED  and  TSEC  is  permissible,  but  for  highest 
precision  it  is  recommended  that  JED  be  the  most  recent  midnight  at  or  before 
the  requested  epoch  and  that  TSEC  be  the  number  of  ephemeris  seconds  elapsed 
between  then  and  the  Interpolation  epoch. 

Output : 


For  I = 1 - 12,  interpolated  quantities  requested  in  IREQ(I)  are  in 
TAB0UT(1,  I)  - TABOUTO,  I)  for  position  vectors  and  in 
T.\B0UT(4,  I)  - TAB0UT(6,  I)  for  velocity  vectors. 

NUT(l)  contains  6iJ;,  the  nutation  in  longitude 

NUT(2)  contains  6e,  the  nutation  in  obliquity 

NUT(3)  contains  6il),  and 

NUT(4)  contains  6e. 

The  units  of  nutations  are  radians  and  radians/time  unit. 

lERR  is  a status  flag. 

lERR  = 0,  normal  return. 

= 1,  epoch  earlier  than  earllst  epoch  on  file. 

= 2,  epoch  later  than  last  epoch  on  file. 

= 3,  IREQ(I)  < 0 or  IREQ(I)  > 2 

“ 4,  Illegal  value  for  ICENT 

Necessary  Modifications:  None 


Reproduced  from  "JPL  Export  Ephemeris  DE-96  Users  Guide",  January  1976. 
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7)  Program  TSTRDE  was  employed  to  exercise  the  binary  ephemeris  file 
using  subroutines  READE  and  SECCAL.  Output  was  matched  to  a microfiche 
copy  of  a JPL  run. 


8)  System  II  of  the  JPL  Ephemeris  Processing  System  produces  an  annual 
tape  for  Prospect  Hill  Radio  Observatory.  This  program  was  modified  to 
operate  with  the  new  JPL96B  binary  ephemeris  file.  Code  for  reading  and 
interpolating  the  obsoleted  binary  file  was  removed  from  subroutines 
EPHEMI,  EPHEl  and  INTRPL  of  System  II  and  replaced  by  one  call  to  READE  in 
EPHEl  along  with  the  required  labelled  common  and  double  precision  statements. 
Provision  for  entering  ET-UT  in  seconds  was  added. 

This  revised  System  II  is  operational  and  has  been  used  satisfactorily  to 
produce  PHRO  ephemeris  tapes  for  1977  and  1978.  Modification  of  System  I 
for  planetary  ephemerides  at  AFGL  can  be  carried  out  when  required. 
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